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Abstract 

We present factorization tlieorems for two exclusive heavy-quarkonium production processes: 
production of two quarkonia in e^e^ annihilation and production of a quarkonium and a light meson 
in S-meson decays. We describe the general proofs of factorization and supplement them with 
explicit one-loop analyses, which illustrate some of the features of the soft-gluon cancellations. We 
find that violations of factorization are generally suppressed relative to the factorized contributions 
by a factor v'^m^Q for each S-wave charmonium and a factor m^Q for each L-wave charmonium 
with L > 0. Here, v is the velocity of the heavy quark or antiquark in the quarkonium rest frame, 
Q = y/s for e~^e~ annihilation, Q = ms for i?-meson decays, is the e^e~ center-of-momentum 
energy, m-c is the charm-quark mass, and is the S-meson mass. There are modifications to the 
suppression factors if quantum-number restrictions apply for the specific process. 
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I. INTRODUCTION 

A crucial step in the calculation of the amplitudes for hard-scattering hadronic processes is 
the separation of the effects of the strong interactions into short-distance and long-distance 
contributions. The short- distance contributions are, by virtue of asymptotic freedom in 
quantum chromodynamics (QCD), perturbatively calculable, while the long-distance con- 
tributions are parametrized in terms of inherently nonperturbative quantities. These sep- 
arations are usually embodied in factorization theorems for the processes. In the case of 
hard-scattering processes that involve heavy-quarkonium states, it has been proposed that 
the effective theory nonrelativistic QCD (NRQCD) could be used to describe the separation 
of perturbative effects that produce a heavy-quark pair from the nonperturbative effects 
that bring about the evolution of the heavy-quark pair into the quarkonium bound state [1]. 
Recently, progress has been made in understanding factorization issues in inclusive heavy- 
quarkonium production processes [2-5]. However, a proof of factorization to all orders in 
QCD perturbation theory is still lacking for inclusive quarkonium production. In the present 
paper we discuss factorization for exclusive quarkonium production. 
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The exclusive production of double-charmonium states in e'^e" annihilation has pro- 
vided an important testing ground in which to compare predictions of theoretical models 
of charmonium production with experimental measurements. Measurements of the cross 
sections for double-charmonium production by the Belle [6] and BABAR [7] collaborations 
have, in several instances, disagreed with theoretical predictions [8-11] and have led to a 
re-examination of the bases for those predictions. 

The exclusive decays of B mesons into a light meson plus a charmonium state are also 
of interest, partly because they could provide new constraints on the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix and enhance our understanding of the origins of CP violation. 
However, the nonperturbative effects of the strong interactions are significant in such pro- 
cesses and must be taken into account in order to make reliable QCD-bascd calculations of 
the process rates. Factorization theorems for these processes would provide a first-principles 
framework within which to take into account the strong-interaction effects. In the case of ex- 
clusive decays of B mesons into a light meson plus a charmonium state, several factorization 
theorems have been proposed [12-14]. 

In this paper, we present proofs, valid to all orders in QCD perturbation theory, of fac- 
torization theorems for the exclusive quarkonium-production processes mentioned above, 
giving details of the proofs that were summarized in Ref. [15]. These are the first proofs of 
factorization theorems for quarkonium production. We also present explicit calculations at 
one-loop order that illustrate key features of the general arguments. Although our analyses 
are for the specific cases of B decays and e'^e~ annihilation, the techniques that we describe 
should apply to other exclusive quarkonium-production processes, and may also shed light 
on factorization in inclusive quarkonium production. However, we note that, because we 
consider exclusive two-body quarkonium-production processes, rather than inclusive quarko- 
nium production, we avoid the issues raised in Ref. [2, 3] concerning light particles that are 
comoving with a quarkonium and the issues raised in Ref. [4, 5] concerning the color-transfer- 
enhancement mechanism that appears when an additional heavy quark is comoving with a 
quarkonium. 

In the analysis of Ref. [15], it was assumed that gluons cannot have transverse momentum 
components that are smaller than the QCD scale, Aqcd- The possibility that external on- 
shell lines can emit gluons of arbitrarily low energy was discussed in detail in Ref. [16], and 
it will be considered here as well. The factorization theorems stated in Ref. [15] remain 
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unchanged. 

In the case of the exclusive production of double-charmonium states in e'^e~ annihilation, 
we will argue that the production amplitude can be written in the following factorized form: 

A (e+e- ^ 7* ^ //i + H^) ^Yl^^^ (^i|Ci|0) (//2|O,|0) . (1) 

The factors (if„|Ci|0) are NRQCD matrix elements, which describe the nonperturbative 
evolution of the charm-quark and the charm-antiquark (cc) pair into a charmonium state 
Hn- The sum over the matrix elements is organized as an expansion in powers of v, the 
relative velocity between the c and the c in the charmonium rest frame. (For charmonium, 

fa 0.3.) The quantity Aij is a short-distance coefficient, which contains the amplitude 
for an e^e~ pair to annihilate through a virtual photon into two cc pairs in the color and 
angular-momentum states of the NRQCD operators Oi and Oj. 

In the case of e"'"e~ annihilation, we define the hard-scattering scale Q = y/s to be the 
center-of-momentum (CM) energy of the e'^e" pair. We will argue that the factorized form 
in Eq. (1) holds up to corrections of relative order /1/2, where // — v^rric/Q for an 5'-wave 
charmonium Hi, fi — rric/Q for an L-wave charmonium Hi with L > 0, and rric is the 
charm-quark mass. As we will discuss in detail, these suppression factors are modified if 
quantum-number restrictions apply for the specific process. 

In the case of exclusive decays of B mesons into a light meson plus a charmonium state, 
we will argue that the decay amplitude can be written in the following factorized form^: 

A{B^H^ + K) = Y,Ff^r^J^{Ml)A,,{H,\0,\Q) + ^ ® ^kj ® ^bi{H,\0,\Q). (2) 

ife ije 

Again, the factors (iJi|(9j|0) are the NRQCD matrix elements, which describe the nonper- 
turbative evolution of the cc pair into a charmonium state Hi. The quantities Fj^^^ , ^kj, 
and $Bi are also nonperturbative objects, which we describe below. The quantities A^^ and 
A'-j^ are short-distance coefficients. They contain the amplitude for the electroweak vertex 
to produce a cc pair in the color and angular- momentum state of Oi. We approximate the 
electroweak vertex as a local four-fermion vertex. The sum over e is over the various oper- 
ators in the electroweak effective action. The sum over / is over the allowed form factors 

^ The two terms in the factorization formula (2) are analogous to the two terms in the factorization formula 
in Eq. (4) of Ref. [12] for the case of decays into two light mesons. 
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that result from shrinking the hard subdiagram (to be described later) to a local vertex 
with respect to the S-meson-to-light-meson transition process. The symbol (8) represents 
the convolution of a short-distance coefficient with the hght-cone distributions of the hght 
meson and the B meson. 

The first term of Eq. (2) contains a 5-meson-to-light-meson form factor 

FfitefiMl) = {mrf^,^)%\B). (3) 

Here, B and K denote the B meson and the light meson, respectively, and Mi is the 
charmonium mass. The quantity I'f{i,e) is the product of a Dirac matrix and a color matrix 
that arises when one shrinks the hard-scattering subdiagram to a point with respect to the 
5-meson-to-light-meson transition amplitude. It is understood that the fields and 
are in a color-singlct state. Following Ref. [12], we define F^^^ in the first term in Eq. (2) 
to be the "physical" meson form factor, which contains both hard and soft contributions. 
Then, in the second term in Eq. (2), one must omit from the short-distance coefficients the 
hard contributions that are already contained in the first term in Eq. (2). 

The second term of Eq. (2) involves the light-cone distribution amplitude(s) of the light 
meson ^Kj, which are defined by the expression 

^ dx+eM-iC^y - l)PKx'-]{K{pK)\^a{x+)P[x+, -x+]*^(-x+)|0) 

J-OC 

= J2^KAy)[^Kjl^, (4) 
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and the light-cone distribution of the B meson $^1, which is defined by the expression 

dx-exp[t^p+x-]{0\^ifs{x-)P[x-,Q]^bamB{pB)) 

X) 

m 

«^ {(A + rribH + fi-<^B2m}ap ■ (5) 

Here, \E' is the quark field, a and /3 are Dirac indices, ^ and ^ in each matrix element 
arc understood to be in a color-singlet state, and the T^j and the Tbui are Dirac-matrix 
structures for the light meson and the B meson, respectively.^ We define light-cone variables 



^ For example, for the leading-twist distributions of the pseudoscalar meson P, the longitudinally polarized 
vector meson V, and the transversely polarized vector meson V±, Tkj is i(/p/4)^75, — ?(/y/4)^, and 
—i{fv±/&)\l^,^], respectively. Here, fp, fv, and fvj_ are the meson decay constants. 
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k — {k'^,k~,k±) in terms of Cartesian components as /c+ = (l/-\/2)(/c° + k^) and k~ — 
(l/-\/2)(A;° — k^). In Eq. (5), n_ is the vector n_ = '\/2(0, 1, 0_l), and we have retained only 
the leading-twist B meson hght-cone distributions. We take the spatial components of pk 

to lie along the minus z direction, and we take the B meson to be at rest. The expression 
[y, x] in Eqs. (4) and (5) is the exponentiated line integral of the gauge field: 

(6) 

P indicates path ordering, is a generator of color SU(3), and is the gluon field. 

In the case of i?-meson decays, we define the hard-scattering scale Q to be the 5-meson 
mass TJiB- We will argue that the factorized form in Eq. (2) holds up to corrections of relative 
order /i, where /i — v^rric/Q for an 5"- wave quarkonium Hi and /i = rric/Q for an L-wave 
quarkonium Hi with L > 0.^ As in the e+e~-annihilation case, these suppression factors 
are modified if quantum-number restrictions apply for the specific process. This result was 
suggested previously in Ref. [12]. However, there it was conjectured only that the violations 
of factorization vanish in the limit rric — t- 0. 

The remainder of this paper is organized as follows. In Sec. II, we specify models for the 
production amplitudes. In Sec. Ill, we outline the proofs of factorization for the processes 
under consideration. There we describe the momentum regions that are leading in the 
hard-scattering scale, the momentum regions in which loop integrands become singular, the 
diagrammatic topologies of the leading and singular regions, the approximations that are 
appropriate to contributions involving momenta that arc soft or collinear, the factorization 
of the soft and collinear singular regions, the subsequent construction of the factorized 
form, and the corrections to the factorized form. We illustrate general features of the 
factorization proof with exphcit one- loop examples in Sec. IV. In Sec. IV A we describe 
the implementation of the soft approximation at the one-loop level. Sections IV B and IV C 
contain one-loop examples for e'^e~ annihilation and B decays, respectively. We summarize 
and discuss our results in Sec. V. The Appendix contains the expressions for the quark- 
antiquark spin-projection operators that we use in Sec. IV. 

^ Ref. [17] presents an analysis of the process B — >■ XcjK in the limit rub — ?■ oo with rric/mb fixed, where nib 
is the bottom-quark mass. The use of the term "factorization" in that paper has, therefore, a different 
meaning than in the present paper, in which we take rric/mb to be a small parameter. 



[y,x]= exp 
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II. MODEL FOR THE AMPLITUDE 



We carry out our analyses in the rest frame of the B meson and in the CM frame of 
the e'^e~ pair, choosing the three-momentum of the quarkonium Hi to be in the positive 
z direction and choosing the three-momentum of the hght meson K or the quarkonium H2 
to be in the negative z direction. We take the constituents of each meson and quarkonium 
to be on the mass shell. We also assume that, for each meson and quarkonium, there is 
an integration over the relative momentum of the constituents, weighted by a meson wave 
function, and subject to the mass-shell constraints. 

We model the B meson as an on-shell "active" bottom quark, which participates in 
the electroweak interaction, and an on-shell "spectator" light antiquark, which does not 
participate in the electroweak interaction. We take the quark and antiquark to be in a 
color-singlet state. We take the bottom quark to have momentum pf, and mass rrib, with (in 
Cartesian coordinates) 

Pb = (^\/ml + ql, Qb^ ~ {mb, Aqcd) • (7) 
We take the spectator antiquark to have momentum pi, with 

Pi^{\Qb\, -Qb) ~ (Aqcd, Aqcd) ■ (8) 

The momentum of the B meson, ps, is given by the sum of p^ and pf. 

Pb ^Pb+Pi = {mB,0) = {\lrnl + q^ + I^bI,©^ ~ (^6,0). (9) 

Similarly, we model the light meson K as an on-shell active light quark and an on-shell 
spectator light antiquark, with the quark and antiquark in a color-singlet state. We can 
write the quark momentum, p^^, and the antiquark momentum, p^^, as 

1 , , 

Pfe, = I^Pk + rfe, (10a) 

Pfc,- = \pk -Tk, (10b) 
with pk ■ Tk — 0. In the rest frame of the light meson, we denote the vectors that are 



8 



associated with the hght meson with a hat. Then, we have 

PK^{mK,0)^{2\n\,0), (11a) 

rk^{0,h), (lib) 

Pk, ^ {\rk\,rk) , (11c) 

Pk, ^ {\rk\,-rk) ■ (lid) 

The quantity is of order Aqcd- 

The boosts from the hght-meson rest frame to the S-meson rest frame are given, for an 
arbitrary momentum k, by 

^^ ^^-^^^ k\ (12a) 
rriK 

k~^^^±^k-, (12b) 
rriK 

k±^ k±. (12c) 
Here, Pcm is the magnitude of the three- momentum of either Hi or K in the 5-rest frame, 

Pcm = 7^ "^6 > (13a) 

Z^y s 

X{x, y, z) = + + 2;^ - 2{xy + yz + zx) , (13b) 
Ejc is the energy of the meson K, 



Ek = y Pcm + mlr^mb. (14) 

Ml ~ rric in Eq. (13a) is the heavy-quarkonium mass, which we will define in our model 
below. Therefore, in the 5-meson rest frame we have 

Pk = [Ek - Pcm] ' ^ + ^cm] , 0^^ ~ (^^^, m^, 0^^ , (15a) 

ffl Ek - Pcm n Ek + Pcm . \ f ^qcB . \ 

, — 7= , Tfcx ~ rrib, Aqcd ■ (15b) 

Vv2 rriK V2 rriK J \ rrib ^ J 

It is now convenient to define a momentum fraction y and a vector q^; that has zero minus 
component. In terms of these quantities, the momenta of the quark and the antiquark are 

Pkg = yPK + qk, (16a) 
Pkq = (1 - y)PK -qk = yPK - Qk- (16b) 
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y is the fraction of minus component of the momentum of the meson that is carried by the 
quark: 



Hence, 



* ^ ( i - .) P. + . (2^ 0, ^ 0, A,eo 1 . (18) 



Finally, we model the charmonium states as an on-shell charm quark and an on-shell 
charm antiquark in a color-singlct state, with the momentum of the charm quark equal to 
Pig and the momentum of the charm antiquark equal to pig. We take 

Piq = ^Pi + Qi, (19a) 
Piq = \Pi - (Ih (19b) 

where Pi is the quarkonium momentum and Pi • qi — 0. In the quarkonium rest frame, we 
denote vectors that are associated with the quarkonium with a hat. The quantity Qi has 
only spatial components, whose magnitudes are of order rricV. Hence, 

Pi = {Mi, 0) = (^2^ml + qf, 0^ , (20a) 
qi = {Q,qi). (20b) 

In the case in which the quarkonium i is in a spin-triplet state, we also define a spin- 
polarization vector ej. In the quarkonium i rest frame, has spatial components of order 
unity and temporal component zero: 

ii = {0,ei), (21) 

which implies that Pj ■ = 0. 

The boost from the rest frame of the quarkonium with momentum Pi to the P-meson 
rest frame or the e~^e~ CM frame is 

k-^^^±^k\ (22a) 

k-^^^^P^k-, (22b) 
k± k±. (22c) 
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The boost from the rest frame of the quarkonium with momentum P2 to the e"^e CM frame 
is 

f + E2- PcM f+ X 

^ ^ M, ^ ' ^^^^^ 

(23b) 

M2 ^ ^ 

^x. (23c) 

Here, 



Ei^^P^^ + Mf^Q, (24a) 

AV^(.,Mf,M|) 
-^^CM = Q , (24b) 



Mi^2^ml + ql (24c) 

M2 = M2 in the case of e+e~ annihilation into two quarkonia, and M2 = ttik in the case of 
i?-meson decays. It then follows that, in the e^e~ CM frame or the 5-meson rest frame. 



p- O ' 



Pl± = 0, P2± = 0, 

ql ~ ^^<5, ~ (25) 

— vmi — ^ 



i m. 



— iit>c — — 

1 O ' ^2 



m, 



qi± ~ ^cV, ei_L ~ 1. 



III. PROOF OF FACTORIZATION 



A. Strategy 

If we dress the lowest-order decay and production amplitudes in our models with addi- 
tional gluons, then certain regions of integration of the gluon momenta yield contributions 
that are leading in powers of the large momentum scale, Q. We will describe these regions 
in Sec. Ill B below. We wish to isolate the contributions from the loop integrations that can 
be calculated in perturbation theory from those that cannot. That is, we wish to isolate 
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contributions in which propagators have large virtuahty, of order Q, from contributions with 
lower virtualities. We call the large- virtuality part of the amphtude the "hard" part. In 
order to establish factorization, we will show that the low-virtuality contributions either 
cancel or can be absorbed into nonperturbative functions. The nonperturbative functions 
are the NRQCD matrix elements for the charmonia and, in the case of 5-meson decays, the 
i?-meson-to-light-meson form factor, the light-cone distribution amplitude for the B meson, 
and the light-cone distribution amplitude for the light meson. We will first demonstrate 
a factorization involving quarkonium distribution amplitudes. Then, we will argue that 
the distribution amplitudes can be straightforwardly decomposed into a sum over NRQCD 
matrix elements multiplied by short-distance coefficients.^ After the factorization of low- 
virtuality contributions, the remaining hard part will depend only on the momenta and 
spins of the quarks and antiquarks that enter into the leading-order process and will be 
independent of the low-virtuality properties of the external mesons. 

The low-virtuality contributions arise from regions of loop integration that are loga- 
rithmically enhanced. In these logarithmically enhanced regions, loop integrations have 
logarithmic power counts and can lead to actual infrared (IR) divergences or would-be IR 
divergences that are cut off by scales smaller than Q, such as quark masses. In the case 
of a would-be divergence that arises from the emission of a gluon that is nearly coUinear 
to one of the external charm quarks, the minimum virtuality of the quark propagator is of 
order m^|fe|/(2|p|), where k is the gluon momentum and p is the charm-quark momentum. 
Hence, the virtuahty can be much less than Q^, and even of order A.qcd- Therefore, we must 
factor such contributions from the hard part in order to arrive at a perturbatively calculable 
contribution. 

One could, in principle, deal with the low-virtuality contributions by devising a suitable 
subtraction scheme for the contributions that would appear order by order in perturbation 
theory. That would be a formidable task, as one would need to ensure that all such contri- 
butions are accounted for in an arbitrarily complicated Feynman diagram, with no double 
counting of contributions. 

For our purposes, we can take a simpler approach. We consider the singularities that 

* For a discussion at the one-loop level of the decomposition of quarkonium light-cone distribution ampli- 
tudes into a sum over NRQCD matrix elements see Refs. [18, 19] . 
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appear in the limit mc — > with qi fixed. First, we estabhsh that the contributions from 
infinitesimal neighborhoods of these singularities can be factored into nonperturbative func- 
tions. Then, we restore rric to its physical value and extend the regions contained in the 
nonperturbative functions from the infinitesimal neighborhoods of the singularities to regions 
of finite size. Then the hard part, which is defined to be the remainder of the amplitude, 
contains no logarithmically enhanced contributions. 

The factorization proofs entail the use of soft and coUinear approximations, which are 
exact at the singular points. These approximations are described in Sees. Ill E and III F. The 
actual factorization is achieved through the use of decoupling relations, which are based on 
the graphical Ward identities of QCD. These decoupling relations are described in Sec. IIIG. 

We note that, because our models make use of on-shell external quarks and antiquarks, 
it is possible to emit coUinear and nearly coUinear gluons of arbitrarily low energy from the 
external lines. This situation is discussed in detail in Ref. [16]. It is unphysical since, in 
a meson, confinement cuts off gluon energies at values of order Aqcd- Nevertheless, it is 
important to establish factorization in the on-shell case in order to guarantee the consistency 
of perturbative calculations of the hard part, which are usually carried out in the context 
of on-shell amplitudes. Because the logarithmically enhanced contributions in the presence 
of a cutoff of order Aqcd are a subset of the logarithmically enhanced contributions in the 
case of on-shell external lines, the factorization argument that we will present also applies 
in the simpler case of a model with a cutoff. As we will see, the methods that we use to 
prove factorization apply to models in which the external particles are off their mass shells, 
provided that the models maintain gauge invariance. For example, one could model the B 
meson as an elementary, color-singlet pseudoscalar that produces the constituent quark and 
antiquark off their mass shells through a pointlike pseudoscalar-interaction vertex that is 
proportional to 75. 

B. Leading momentum regions 

In describing the regions of loop momenta that yield contributions that are leading in 

powers of the large scale Q, we make use of the nomenclature of Ref. [16]. We first describe 
the various regions of momentum space, and then, in Sec. IIIB5, we specify the conditions 
that must be fulfilled in order for these regions to give leading contributions to an amplitude. 
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In a Feynman diagram, we call a gluon or quark, or, generically, a line that carries momentum 
of type X an "X gluon," "X quark," or "X line." 

1. Hard, soft, and collinear regions 

The hard {H), soft (S), collinear-to-plus (C^), and coUinear-to-minus (C~) momenta 
have components with the following orders of magnitude: 

H: g(l,l,U), (26a) 

S: ge5(l,l,l±), (26b) 

C+: ge+[l,(^+)^r7+], (26c) 

C-: g6-[(7;-)^l,r/I]. (26d) 

The energy scales of the various types of momenta arc determined by the parameters e^, 
e"*", and e~. The soft region of momentum space is defined by the condition 

es « 1. (27) 

The collinear regions of momentum space are defined by the conditions 

e± < 1 

< 1. (28) 

In the case of S-meson decays, there is also a leading region that is associated with 
momenta that are nearly collinear to the light-quark momentum pi. We call this region CK 
It is characterized by momenta that scale as 

C': Qe% + iv'fei + v'en, (29) 

where e^, is a unit vector that is parallel to the lightlike vector pi, ei is the parity inverse of 
ei, and e[ is a unit vector that is transverse to e; and q. The C' region is defined by 



e 



I < 



Aqcd/ Q, 



rf < 1. (30) 
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We assume that pi does not lie exactly in the plus or minus direction. In order to simplify 
the discussion to follow, we often do not mention C' momenta explicitly. In these instances, 
it may be assumed that the lines carrying C' momenta may be treated analogously to the 
lines carrying momenta. 

Wc note that soft and collinear contributions lie in restricted regions of phase space. 
Were it not for enhancements that arise from propagators with low virtuality, soft contribu- 
tions would be suppressed by a phase-space factor and collinear contributions would be 
suppressed by a phase-space factor (e='=)^(77^)^. The low-virtuality propagators associated 
with these contributions lead to loop integrals that have a logarithmic power count and to 
contributions that behave as and (e^)''(?7^)°. We refer to such contributions as soft and 
collinear logarithmic enhancements. 

The definitions given above for the H, S, and C* momentum regions do not specify 
unambiguously the boundaries between them. For instance, if the 77* parameters in the 
collinear regions take on values that are not too different from one, then the momenta 
are not distinguished from the S momenta; i.e., it would not be clear if a momentum with 
77' close to 1 belongs to the C* region or to the S region. Hence, the possibility of double 
counting arises. Analogous issues appear at the other boundaries between the H, S, and C* 
regions. However, as we have mentioned previously, this is not a problem for our proof of 
factorization, which will be presented later, because the proof focuses on the singularities 
and would-be singularities, rather than on the momentum regions. The heuristic description 
of regions presented here is intended only to set the stage for the subsequent discussions of 
the singular regions. 

In contrast with the corresponding momentum regions, the soft and collinear singularities 
are distinct. The soft singularities appear in the limit — )■ or e* — )■ and the collinear 
singularities appear in the limits 77* — > 0. A double (soft and collinear) singularity can arise 
if — > and 77* — > at the same time. 

2. Endpoint region 

In the case of 5-meson decays, there is a leading contribution from the so-called "end- 
point" region [12]. This contribution is associated with a gluon that connects the i?-meson 
and light-meson antiquarks to the remainder of the amplitude. The contribution in the 
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endpoint region arises from a would-be infrared divergence that corresponds to the singular 
point at y = 0. The divergence is cut off by g^, the residual momentum of light- meson 
antiquark, and by pi, the momentum of the meson antiquark, both of which are of or- 
der Aqcd- That is, the divergence is cut off at y ~ Aqcd/'"^6- In our model, the explicit 
diagrammatic factors yield a linearly divergent power count, but the would-be divergence 
is moderated by a factor y from the light-meson wave function and is actually logarithmic. 
The gluon that is associated with the endpoint region carries S momentum of order Aqcd-^ 
If the gluon that is associated with the endpoint region attaches to an active-quark line 
from the B meson or the light meson or to a heavy-quark or heavy- antiquark line from the 
charmonium, then its momentum causes the propagators of those lines to be off shell by an 
amount of order rrift Aqcd- We call such lines "semihard" lines. Contributions from these 
lines can be calculated in perturbation theory. We treat the semihard region as part of the 
hard region, and we include lines carrying semihard momenta in the hard subdiagram that 
we will describe below. 

3. Glauber region 

The "Glauber" region is also leading in power counting [1, 24]. In this region, <^ 
\k~\ <^ and k"^ <^ Q^. In processes with two incoming hadrons, such as Drcll-Yan 
lepton-pair production, pinch singularities can develop in the Glauber region for the and 
k~ contours of integration on a diagram- by-diagram basis [24-26]. The pinches arise when a 
gluon connects a spectator parton in one initial-state hadron with a spectator parton in the 
other initial-state hadron. (Here, in contrast with the terminology that is used to discuss 
exclusive i?-meson decays, "spectator parton" means a parton that does not participate 

^ If the spectator antiquark line that connects the B meson to the hght meson carries a C~ momentum 
whose invariant square is of order Aq(^j-,/Q, then that momentum is said to be in the soft-coUinear or 
messenger region [20] . Such a momentum arises from a small part of the phase space in which a gluon on 
the S-mcson side of the soft-collincar spectator line carries away most of pi and a gluon on the light- meson 
side of the soft-coUinear spectator line carries away most of y»/c- . It has been argued that the soft-coUinear 
region is leading only when one makes use of certain infrared regulators [21-23]. In any case, a contribution 
from the soft-coUinear region does not require any special treatment in our factorization argument: The 
gluon on the light- meson side of the soft-collinear spectator line can be treated as C+, and the gluon on 
the B-meson side of the soft-collinear spectator line can be treated as S, as it would be in the endpoint 
region. 



16 



in the hard-scattering process.) The pinches appear because the momentum of a gluon 
that attaches to a spectator-parton hne must route through the hadron wave function and 
the active-parton hne from that hadron to the hard process. If the gluon's momentum 
in the active-parton hne is in the same direction as the momentum of the active parton, 
then, in the spectator-parton hne, it is in the direction opposite to the momentum of the 
spectator parton. Consequently, there is a pinch in the hght-cone variable that is conjugate 
to the direction of the momentum of the hadron. In contrast, in exclusive processes, all 
of the partons in a hadron are connected in the lowest-order process, either through the 
hard subprocess or, possibly through a soft gluon in the case of S-meson decays. (See 
the discussion of the cndpoint region above.) Thus, if an additional gluon carrying soft 
momentum attaches to a parton, one can always route that momentum through a leading- 
order connection to the hard part, avoiding routings through other partons in the hadron 
that could produce a pinch. Because of this, the k'^ and k~ contours of integration are not 
pinched in the Glauber region in exclusive processes, and it is possible to deform them out 
of the Glauber region on a diagram-by-diagram basis. Therefore, we ignore the Glauber 
region in the remainder of our discussion. 

4- Threshold region 

In the case of a quarkonium, there are "threshold enhancements" that are associated with 
the exchange of a gluon between the quark and the antiquark. (See Ref. [1] for examples.) In 
the quarkonium rest frame, the enhancement occurs when the exchanged gluon has momen- 
tum components ~ rricV^ and |fc| ~ rricV. The enhancement produces a power infrared 
divergence that is cut off by the relative momentum of the quark and antiquark q ~ rricV. 
The divergence is proportional to Tnc/\q\ ~ 1/v. Now let us consider the momentum of 
the exchanged gluon in the e'^e~ CM frame in the case of e'^e~ annihilation and in the 
S-meson rest frame in the case of S-meson decays. In these frames, as can be seen from the 
boosts in Eqs. (22) and (23), an exchanged gluon in the quarkonium with momentum Pi 
has momentum k ~ {Qv,mlv/Q,mcV±), and an exchanged gluon in the quarkonium with 
momentum P2 has momentum k ~ {mlv / Q , Qv , rricV ±) . Therefore, the exchanged gluons 
associated with threshold enhancement have C+ or C~ momentum, and, in our analysis, 
we do not distinguish them from other gluons with C"*" or C~ momentum. Because the 
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threshold enhancements involve the gluons and heavy quarks in a single quarkonium, in 
each Feynman diagram they are a priori compatible with the factorized forms. Therefore, 
it will not be necessary to manipulate the threshold contributions or to identify them by 
considering the limit v ^ 0. 

5. Leading momentum configurations in Feynman diagrams 

Next we identify the configurations of momentum types that can yield leading contribu- 
tions in the Feynman diagrams. By "leading" , we mean contributions that are not suppressed 
as powers of ratios of momentum components. We follow the analysis presented in Ref. [16]. 
Here, and throughout this paper, we work in the Feynman gauge. 

We start with a basic diagram that is just the amplitude of lowest order that involves 
the external quark and antiquark from each meson. Then we add gluons, one at a time, 
determining for each gluon the momentum types that produce leading contributions. The 
added gluons can contain quark, gluon, and ghost vacuum-polarization loops. 

Because there are many redundant ways to obtain a given momentum configuration in a 
diagram, it is useful to define a convention for the way in which we add gluons. In order to 
do that, we first define combination momenta and CC. A momentum arises from 
the sum of a momentum and an S momentum with €5 ~ e^rj^ or from the sum of a C"^ 
momentum and a momenta with e^{r]'^y' <^ <^ e^. A CC momentum arises from the 
sum of a momentum and a momentum with e''" ~ e~. These combination momenta 
have the following orders of magnitude: 



C+r Qe+(l,77+,r7l), (31a) 

C-: Qe-{r,l,Vl), (31b) 

CC: Qecc{l,l,ricc±), (31c) 

where 

1 > > iv^f. (32) 

Analogous combination momenta may be defined for combinations of momenta with other 
momenta. A C' momentum has dominant component in the e; direction. A CC' momentum. 
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ep(?7p ) 


~ fc ~ p p 'P ~ fc ~ ^p 


ecCfc ~ ecCp 



TABLE I: Conditions that a gluon with momentum k must fulfill in order to attach to a line with 
momentum p. These conditions guarantee that the resulting attachment is allowed according to 
the convention described in the text and that it results in a leading contribution. Here, "leading" 
means that the contribution is not suppressed as powers of ratios of momentum components. In 
each table, the left-hand column gives the momentum type of the gluon with momentum k, and 
the top row gives the momentum type of the line with momentum p. The symbol "~" means that 
quantities are of the same order. For purposes of power counting, an H line behaves as a soft line 
with es ~ 1. The rules for attachment when k is C^ are the same as the rules for attachment when 
k is C^. If k is S, and the lines to which it attaches have momentum pi and pj, then pi and pj 
cannot both be C"^ or C~ . If A; is C^, then at least one of pj and pj is C^. Analogous conditions 
exist for the attachments of gluons with C' momenta. 

which arises from the sum of a momentum and a momentum with e' ~ e^, has at 
least two components of order e'. 

Now we define our convention for adding gluons to the basic diagram. We say that a 
gluon with momentum / can attach to a line with momentum p only if the energy scale of 
the momentum p + / is of the same order as the e parameter of the momentum p and one of 
the following conditions is fulfilled: 

1. The momentum p + / is of the same type as the momentum p; 

2. The momentum p is and p + I is C^; 

3. The momentum p is or C" and p + l is CC. 

The analysis in Ref. [16] shows that, if we consider only the terms 2p ■ I in propagator 
denominators, then the gluons that we add to the basic diagram must be S, C'^, C~ , or 
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in order to obtain a leading contribution. The combination momenta defined above, 
arise when we add S, C^, and C' momenta. If we consider, as well, the terms and in 
propagator denominators, then contributions are subleading unless 

k-p>p\ (33) 

The constraints in Eq. (33) lead to additional restrictions on the momentum combinations 
that yield leading contributions. These restrictions, combined with our conventions for 
adding gluons to a diagram, result in the rules for the attachments that yield leading contri- 
butions that are given in Table I. The rules in Table I also apply when the gluon attaches 
to one of the fermion lines that begins as an external quark or antiquark. In that case, one 
sets 77=*= = for the external quark or antiquark. We have not displayed the rules for the 
attachments of gluons with or momenta to lines with or momenta because 
the rules for such attachments are complicated and cannot be characterized simply in terms 
of the magnitudes of the momentum components. For our purposes, it suffices to note that 
necessary conditions for such attachments are given in Eq. (33). 

Some of the allowed attachments in Table I change the type of the momentum in the top 
row, for example, when we add an S gluon to a gluon with eg ~ 77^6^. That change can 
propagate through the Feynman diagram. In those cases one must check that the rules in 
Table I still allow the attachments of all the vertices that are affected by the change. 

The constraints in Eq. (33) imply that an attachment of a gluon to a given line is allowed 
only if the virtuality that it produces on that line is of order or greater than the virtuality 
that is produced by the gluons that attach to that line to the outside of the attachment in 
question. If a gluon with momentum k of type C^, C^, S, C^, or CC attaches to a hue 
from an on-shell external quark or antiquark, then it adds virtuality Q'^^kiVik)'^ i Q'^^k^ki 
Q'^es^, Q^Cfc , or Q^ecCfe, respectively. 

C. Topologies of the leading regions 

Now let us specify the diagrammatic topology that corresponds to the leading regions. 
The topology of the leading regions for e'^e~ annihilation into two quarkonia is shown in 
Fig. 1. In this topology, there is a hard subdiagram that includes the lowest-order process, a 
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FIG. 1: Leading regions for double-charmonium production in e^e annihilation. The wavy hne 
represents the virtual photon. 

soft subdiagram, and a jet subdiagram for each of the two collinear regions, which correspond 
to the two quarkonia. In the hard subdiagram, all propagator denominators are of order 
s. The soft subdiagram includes gluons with soft momenta and loops involving quarks and 
ghosts with soft momenta. The soft subdiagram attaches to the jet subdiagrams through 
any number of soft-gluon lines. Each jet subdiagram contains the quark and antiquark 
lines for a given quarkonium, as well as gluons and loops involving quarks and ghosts with 
momenta collinear to the meson or quarkonium. The subdiagram attaches to the hard 
subdiagram through the quark and antiquark lines and through any number of gluons. 
As was pointed out in Ref. [16], because gluons with momenta of arbitrarily low energy 
can contribute at leading power in Q, the subdiagram also attaches to the soft and 
subdiagrams through any number of gluons. 

There are two distinct topologies in the case of 5-meson decays: one in which the B- 
meson and light-meson spectators participate in the hard interaction and another in which 
they do not. These two topologies are shown in Figs. 2(a) and 2(b), respectively. The 
topology of Fig. 2(a) is appropriate when the light-meson- antiquark momentum is outside 
the endpoint region, and the topology of Fig. 2(b) is appropriate when the light-meson- 
antiquark momentum is in the endpoint region. 

In each topology in Fig. 2, there is a hard subdiagram that includes the lowest-order 
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parton-level process, there is a soft subdiagram, and there is a jet subdiagram for each of 
the two coUinear regions, which correspond to the hght meson and the quarkonium. In 
the hard subdiagram, all propagator denominators are of order ml or mftAgcD- The soft 
subdiagram includes gluons with soft momenta and loops involving quarks and ghosts with 
soft momenta. The soft subdiagram attaches to the jet subdiagrams and to the i?-meson 
quark and antiquark quark lines through any number of soft gluon lines. The J"*" subdiagram 
contains the quarkonium quark and antiquark lines; the J~ subdiagram contains the light- 
meson quark and antiquark lines; the J' subdiagram contains the S-meson spectator-quark 
line. In addition, the jet subdiagrams contain gluons and loops involving quarks, gluons, and 
ghosts with momenta in the C^, or regions. Each jet subdiagram contains the active- and 
spectator-quark lines for a given meson or quarkonium as well as gluons and loops involving 
quarks, gluons, and ghosts with momenta coUinear to the meson or quarkonium. A J"^ or J' 
subdiagram attaches to the hard subdiagram through the active- and spectator-quark lines 
in the topology of Fig. 2(a), through the active-quark lines in the topology of Fig. 2(b) and 
through any number of gluons. (We have not shown explicitly the attachments of the J' jet 
subdiagram that involve gluons with C' momenta.) As we have already mentioned, because 
gluons with C* momenta of arbitrarily low energy can contribute at leading power in Q, the 
subdiagram also attaches to the soft, J^, and subdiagrams through any number of 
gluons, and the subdiagram also attaches to the soft and subdiagrams through 
any number of gluons [16]. 

In the case of the topology of Fig. 2(b), we show exphcitly a gluon that is marked with 
an asterisk. This is the gluon that was mentioned in our discussion of the endpoint region 
in Sec. IIIB2. We choose the momentum routing so that it always carries the momentum 
of the 5-meson antiquark and the (endpoint) momentum of the light-meson antiquark, 
both of which are of order Aqcd- Therefore, we consider this gluon to be part of the 
soft subdiagram. However, we single it out because it must be present in our model in 
order for the light antiquarks (spectators) to be connected to the remainder of the diagram 
and because its momentum is fixed by the S-meson and light-meson antiquark momenta. 
Soft gluons and low-energy gluons can connect to the marked gluon, although we have 
not shown these connections explicitly. We show the marked gluon connecting to the hard 
subdiagram because its allowed connections to the jet subdiagrams or the 6-quark line result 
in propagators with semihard virtualities, of order mbAqcB, which are part of the hard 
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(a) (b) 

FIG. 2: Leading regions for the S-meson-decay case. The coUinear region J"*" corresponds to the 
charmonium and colhnear region J~ corresponds to the hght meson. 

subdiagram. As we have said, the topology of Fig. 2(b) appUes to the endpoint region, in 
which the marked gluon has virtuaUty of order Aq^q. Away from the endpoint region, the 
marked gluon itself has virtuality of order mftAqcD and can be incorporated into the hard 
subdiagram, resulting in the topology of Fig. 2(a). 

D. Topologies of the Singular Regions 

In the massless sector of QCD, there are singularities that are associated with soft and 
coUinear divergences. (See, for example, Refs. [25-28].) In the present case, the masses of 
charm quarks and antiquarks cut off some of the colhnear divergences. Some potential soft 
divergences are also cut off because they are associated with a gluon that attaches to a line 
that cannot go precisely to its mass shell because a quark mass cuts off a colhnear diver- 
gence. Nevertheless, as we have mentioned, we wish to consider not only actual divergences, 
but also divergences that appear only in the limit m^-jQ — > with gj fixed, because such 
divergences are associated with logarithmic enhancements. These divergences are associated 
with singularities in the domain of integrations. We call the infinitesimal neighborhoods of 
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such singularities "singular regions." In the remainder of the discussion of the factoriza- 
tion of the contributions of the singular regions, we assume that we have taken the limit 
i^dQ with fixed. 

The topologies of the singular regions follow from the power-counting rules that are given 
in Sec. IIIB. These topologies have been discussed in Ref. [16]. Here, we recapitulate 
that discussion, describing the relationships of the topologies of the singular regions to the 
topologies of the leading regions in Figs. 1, 2(a), and 2(b). 

The singular region is situated in the outermost part of the J* subdiagram. (Here, 
"out" means toward the external fermion lines.) We call this part of the J* subdiagram the 
J* subdiagram. (Wc denote the part of the J* subdiagram that excludes the J* subdiagram 
as the J* — J* subdiagram.) The S singular region is situated in the outermost part of the S 
subdiagram. We call the S singular part of the S subdiagram the S subdiagram. (We denote 
the part of the S subdiagram that excludes the S subdiagram as the S — S subdiagram.) S 
singular gluons connect the S subdiagram only to the J' subdiagrams and to the external 
6-quark line. The gluon that is marked with an asterisk in the topology of Fig. 2(b) is 
not part of the S subdiagram because its momentum components are fixed to be of order 
^QCD- That is, it is S but not S singular. The J* subdiagrams connect to the J*, and 
H subdiagrams via C" gluons. We denote by H the union of all of the subdiagrams in our 
topology except for S, J"*", J~, and J'. We note that the connections of the subdiagrams 
to the S subdiagram via gluons were not considered in the discussions in Refs. [25, 26]. 
Otherwise, the general structure of the topologies of the singular regions that we consider 
are the same as in Refs. [25, 26], provided that we identify the hard subdiagram in those 
references with H. 

E. CoUinear approximation 

We now describe the coUinear approximations, which are useful in factoring the sin- 
gular contributions. We follow the notation of Ref. [16]. 

Suppose that there is a gluon with momentum in the singular region that attaches 
to a line that is not in J*. Then, we can apply a coUinear approximation to that gluon 
[24-26] without loss of accuracy. The C* approximation consists of replacing g^^ in the 
gluon-propagator numerator as follows: 
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Kni, (C+), 



k ■ fix — ie 

{C-), (34) 



k • n2 + is 
kfjjiii, 



^ k ■ hi — is 

Here, the index // corresponds to the attachment of the gluon to the hne with momentum 
not in the singular region, and the index v corresponds to the attachment of the gluon 
to the J' subdiagram. We always use the convention that k flows out of a (7+ or C' hne 
and into a C~ line. There is a large amount of freedom in choosing the auxiliary vectors 
hi, 77,2 and hi in Eq. (34). We need only have hi ■ pig > (or hi ■ piq > 0), n2 ■ P2q > (or 
^2 • P2q > 0) , and hi ■ pi > in order to reproduce the amplitude in the coUinear singular 
region. Our choice is to take hi, fi2, and hi to be lightlike vectors in the minus, plus, and 
minus directions, respectively: 

fii = rii = (l/V2)(0,l,Ox), (35a) 
n2 = ri2 = (l/V2)(l,0,0^), (35b) 
hi = hi. (35c) 

In order for the C* approximation to be exact in the C* limit, j ■ k must be equal to j^k"^, 
where j is the current to which the index of the gluon with momentum k attaches and rij 
is a unit lightlike vector in the direction. This requirement is met provided that the gluon 
does not attach with its /i index to a line that is also carrying momentum in the singular 
region. That is, the C* approximation holds in the coUinear limit if j is a current in any of 
the subdiagrams except for the subdiagram. We note that, in the C" approximation, the 
gluon's polarization is longitudinal, i.e., proportional to the gluon's momentum. This fact is 
essential to the apphcation of graphical Ward identities to derive decoupling relations. We 
note also that the coUinear approximation is exact, not only for the coUinear singularity, but 
also for the associated coUinear logarithmic enhancement. 



F. Soft approximation 



We now describe the soft approximation, which is useful in factoring the S singular 
contributions. Again, we follow the notation of Ref . [16] . 
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Suppose that there is a gluon with momentum k in the S singular region that attaches to 
a hne carrying momentum p that hes outside the S singular region. Then we can apply the 
soft approximation to that gluon without loss of accuracy. The soft approximation [29, 30] 
consists of replacing g^y in the gluon-propagator numerator as follows: 



where the index jj, corresponds to the attachment of the gluon to the line with momentum 
P- 

Unlike the collinear approximation, the soft approximation depends on the momentum 
of the line to which the gluon attaches. However, it is convenient to apply the same soft 



collinear either to the momentum of the quark or the momentum of the antiquark in the 
jet. In the case of the light-quark jet, the quark and antiquark momenta are parallel, up to 
corrections of relative order Aqcd/Q- In the case of the quarkonium jet(s), the quark and 
antiquark momenta Piq and Piq are parallel up to corrections of relative order rricv/Q. In 
both cases, we neglect the difference between the quark and antiquark momenta and define a 
"modified soft approximation" for each jet that corresponds to the soft approximation for the 
average of the quark and antiquark momenta. The leading errors that arise in applying the 
modified soft approximation to the light- meson, charmonium-1, and charmonium-2 jets are 
of relative order Qk^/p^ ~ Aqcd/<5, 9i±/-Pi" ~ rricv/Q, and q2A_/ P2 ~ rricv/Q, respectively. 

It is convenient, for purposes of discussing the decoupling relations in Sec. Ill G, to choose 
lightlike vectors for the soft approximation that correspond to the average of the quark and 
antiquark momenta in the limits Aqcd/Q ~^ ^^^nd nic/Q 0. In making this choice, 
we introduce an error of relative order Aqq^/Q^ in the case of the light-meson jet and of 
relative order m^JQ^ in the case of the quarkonium jet(s). These errors are neghgible in 
comparison with the errors that we make in neglecting the difference between the quark and 
antiquark momenta. Then, for both the light-meson and quarkonium jets we have the same 
soft approximations. 

For the attachment of the gluon with momentum k to any line with momentum in the 
C"*" (C~) singular region, the (modified) soft approximation consists of the following replace- 




■P 



(36) 



approximation to all of the lines in the subdiagram. The lines in the subdiagram are 
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ments in the gluon-propagator numerator: 



(37a) 



k ■ rii + i£ 



is-), 



(37b) 



k ■ n2 — is 



where rii is a hghthke vector that is proportional to Pi and n2 is a hghthke vector that is 
proportional to P2 or px- We normalize rii and n2 so that they are the parity inverses of 
the vectors ni and n2 in Eq. (35), respectively: 



The index n contracts into the line carrying the momentum of type C"^ (C*~)- As we have 
mentioned, the modified soft approximation in Eq. (37) accounts for the contributions in the 
S singular region up to corrections of relative order Aqqb/Q in the case of the light-meson 
jet and up to corrections of relative order rricv/Q in the case of quarkonium jets. We note 
also that the soft approximation is valid at this accuracy not only for the soft singularity, 
but also for the associated soft logarithmic enhancement. 

We do not apply soft approximations to the B meson because the eikonal vectors that 
are associated with soft approximations for J' and the 6-quark line are not approximately 
proportional to each other. That is, a common soft approximation cannot be applied to 
the B meson. In consequence, the cancellations of the soft contributions that apply in the 
case of the quarkonia and the light meson (described in Sec. Ill 1 1) fail in the case of the B 
meson. 

G. Decoupling relations 

Once we have implemented a coUinear approximation or a soft approximation, the associ- 
ated gluons are longitudinally polarized. This allows us to make use of decoupling relations 
to factor gluons with momenta in the soft or coUinear singular regions from certain parts of 
the amplitude. The general graphical form of the decoupling relations for longitudinally po- 
larized gluons is shown in Fig. 3. A decoupling relation of this form applies when any number 
of longitudinally polarized gluons attach to a subdiagram in all possible ways, provided that 





712 = (1/V2)(0,l,0x). 



(38b) 
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FIG. 3: Graphical representation of the decoupUng relations for collinear gluons and the decoupling 
relations for soft gluons. The applicability of these decoupling relations is described in the text. 
The relations show the decoupling of longitudinally polarized gluons, which are represented by 
curly lines. The longitudinally polarized gluon lines are to be attached in all possible ways to 
the Green's function that is represented by an oval. The factors k^n'^ /{k ■ fn) \k^n^/{k ■ rii)] that 
appear in the collinear (soft) approximations are represented by the arrows on the gluon lines. 
The external lines with hash marks are truncated. In addition, the subdiagram can include any 
number of untruncated on-shell external legs, provided that the polarizations of the on-shell gluons 
are orthogonal to their momentum. The pi are momenta, and the ai are color indices. The double 
lines are C~^, C~ , 5+, or S~ eikonal lines, as is described in the text. 

the gluon momenta are all proportional to each other. ^ 

If the external gluons all have momenta in one of the C* singular regions, then the gluon 
momenta are all proportional to each other, and a decoupling relation of the form in Fig. 3 
holds, once the C* approximation has been implemented to render the gluon polarizations 
longitudinal. The subdiagram can have any number of truncated legs and any number of 
untruncated on-shell external legs (not shown in the figure) , provided that the polarizations 
of the untruncated on-shell gluons are orthogonal to their momentum. The eikonal (double) 
lines in this decoupling relation have the Feynman rules in the C~^, C~ , or C' cases that a 
vertex is ^igTafiif^, ±igTan2fj., or igTafiifj^ and a propagator is i/{k ■ fii — is), i/{k ■ n2 + ie), 
i/{k-ni — ie), respectively, where the upper (lower) sign in the vertex is for eikonal lines that 
attach to quark (antiquark) lines. Here, is an SU{3) color matrix in the fundamental 
representation. (Our convention is that a QCD gluon-quark vertex is igTa'Ju.) We call these 
eikonal lines "C* eikonal lines." The Feynman rules for the eikonal lines in these decoupling 

^ In the case of an Abelian theory, such as quantum electrodynamics, a decoupling relation of this form 
holds even if the gluon momenta are not proportional to each other. 
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relations are summarized in the first, second, and third hnes of Table II. 

If the external gluons all have momenta in the S singular region, then the decoupling re- 
quirement that the gluon momenta be proportional to each other is not necessarily satisfied. 
However, if the subdiagram into which the S singular gluons enter is a {J~) subdiagram, 
then a soft momentum k entering that subdiagram contracts only into currents proportional 
to rii (712), up to corrections of relative order qk2/PK ~ ^qcb/Q for a light- meson jet and 
relative order qi±/ ~ Q2±/ P2 ~ rricv/Q for a charmonium jet. Consequently, at these 
levels of accuracy, we can make the replacement 

k ^ ki — fii ^^ ^ (39a) 

TT-i • TT-i 

in the J"*" subdiagram and associated soft approximation and the replacement 

k^h^Ui ^^ (39b) 

77,2 ■ n2 

in the J~ subdiagram and associated soft approximation [25, 26]. Since 

rii ■ ki — rii ■ k, (40a) 
n2-k2 = n2- k, (40b) 

these replacements do not change the amplitude, up to corrections of relative order rndQ 
and Aqcd/Q- In subsequent discussions, we consider these replacements to be part of the 
modified soft approximation. After these replacements have been made, the gluon momenta 
entering the J"*" ( J~) subdiagram are all proportional to each other, and a decoupling relation 
of the form in Fig. 3 holds. In these decoupling relations, the eikonal lines have the Feynman 
rules that a vertex is -^igTani^ {^igTan2ij,) and a propagator is i/{k-ni + ie) [i/{k-n2— is)] 
when the subdiagram is C+ (C~). These rules follow from Eq. (40). We call these eikonal 
lines S'^ and S~ eikonal lines, respectively. The Feynman rules for the eikonal lines in the 
soft decoupling relations are summarized in the fourth and fifth lines of Table II, respectively. 

H. Factorization of the singular regions 

Now we summarize the factorization of the singular regions. We refer the reader to 
Ref. [16] for detailed arguments. 
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Type 


Vertex 


Propagator 


C+ 




i 

k ■ fii — ie 


c- 


±igTan2^ 


i 

k ■ n2 + ie 


& 




i 

k ■ ni — ie 


s+ 


±igTani^ 


i 

k ■ ni+ie 


s- 




i 

k ■ n2 — ie 



TABLE II: Feynman rules for the coUinear (C^ and C') and soft (S^) eikonal lines. The upper 
(lower) sign is for the eikonal line that attaches to a quark (antiquark) line. 

In analyzing the singular regions, we wish to identify the momentum configurations that 
yield singular contributions. We can do so by making use of the power-counting rules that 
we have outlined in Sec. IIIB and invoking the following specific interpretations of those 
rules: the symbol ~ and the phrase "of the same order" mean that quantities differ by 
a finite factor, while the phrases "much less than" and "much greater than" mean that 
quantities differ by an infinite factor. It follows that, for gluons in the singular regions, our 
convention that an allowed attachment of a gluon cannot change the essential nature of the 
momentum of the line to which it attaches has the following meaning: The attaching gluon 
cannot have an energy that is greater by an infinite factor than the energy of the line to 
which it attaches. 

The rules in Sec. IIIB lead to complicated relationships between the allowed momenta 
of gluons in a given diagrammatic topology. However, there is a general principle, which we 
have already mentioned, that allows us to organize the discussion: The attachments of gluons 
to a given line must be ordered so that a given attachment produces a virtuality along the 
line that is of order or greater than the virtualities that are produced by the attachments 
that lie to the outside of it. In particular, the virtuality that a C*, or S* singular gluon 
produces on a line with j ^ i or an S line is of order the energy of gluon times the energy 
of the line to which it attaches. 

Our goal is to factor contributions from all subdiagrams except J^, to factor all 
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singular contributions from all subdiagrams except J' and the external 6-quark line and to 
factor all S singular contributions from the subdiagrams. We will show that the factored 
soft contributions that are associated with the external-quark and external- antiquark lines 
in ultimately cancel. 

Note that we do not factor S singular contributions from or from the external b- 
quark line. Nor do we factor singular contributions from the external 6-quark line in 
the S-meson subdiagram. In principle, we could carry out such factorizations. However, 
because the soft approximations are different for the b quark and the light antiquark in the B 
meson, we do not expect the factored soft contributions that are associated with the 6-quark 
and light-antiquark lines to cancel. Furthermore, it will prove convenient, for purposes of 
expressing our results in terms of a i?-meson light-cone distribution, not to factor the 
singular contributions from the 6-quark line in the S-meson subdiagram. 

In the singular limits eg — > 0, 77' — > 0, — > 0, an infinite hierarchy of energy scales 
emerges. The energy scales of the various levels in the hierarchy are separated by infinite 
factors. We characterize each level in the hierarchy by the energy scale of the S singular 
gluons in that level. We call this scale the nominal energy scale of that level. Collinear 
singular gluons in a level may have energies that are of the nominal energy scale or energies 
that are infinitely larger than the nominal scale, but still infinitesimal in comparison with the 
nominal energy scale of the next higher level. We call the latter gluons "large-scale coUinear 
singular gluons." We carry out the factorization iteratively, starting with the level with the 
largest nominal energy scale. As we shall see, this ordering of the factorization procedure 
is convenient because it allows us to apply the decoupling relations rather straightforwardly 
to decouple gluons whose connections lie toward the inside of the Feynman diagrams before 
we decouple gluons whose connections lie to the outside of the Feynman diagrams. 

We will illustrate the factorization of the large-scale collinear gluons and the nominal-scale 
soft and collinear gluons for the case of double-charmonium production in e'^e" annihilation 
by referring to the diagram that is shown in Fig. 4. In this diagram, we have suppressed 
gluons with energies that are much less than the nominal scale. These gluons have connec- 
tions that lie to the outside of the connections of the gluons that are shown explicitly. In the 
diagram in Fig. 4, each gluon represents any finite number of gluons, including zero gluons. 
For clarity, we have suppressed the antiquark lines in each meson and we have shown explic- 
itly only the connections of the gluons to the quark line in each meson and only a particular 
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FIG. 4: Diagram to illustrate the factorization of large-scale collinear gluons and nominal-scale 
soft and collinear gluons for the case of double-charmonium production in e~^e~ annihilation. C^g 
denotes a large-scale C* singular gluon, C^g denotes a nominal-scale singular gluon, and Sns 
denotes a nominal-scale S singular gluon. 

ordering of those connections. However, we take the diagram in Fig. 4 to represent a sum of 
many diagrams, which include all of the connections that we specify in the arguments below 
of the singular gluons to the quark and antiquark in each meson, to other singular gluons, 
and to the H sub diagram. 

1. Factorization of the large-scale C* singular gluons 

First, we factor the large-scale C* singular gluons. In the first step of the iteration, these 
include gluons with finite energies, as well as infinitesimal energies. In subsequent steps, only 
gluons with infinitesimal energies are involved. There is a hierarchy in the energy scales of 
the large-scale C* singular gluons. We factor these gluons iteratively, beginning with the 
largest energy scale. 

We apply the C* approximations and the decoupling relations. In applying the 
decoupling relations, we include the attachments that are allowed by our conventions to all 
subdiagrams outside of J^, and, in applying the decoupling relations, we include the 
attachments that are allowed by our conventions to all subdiagrams outside of and the 
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external &-quark line. We also include, formally some attachments that may yield vanishing 
contributions in the singular limits. These are attachments to H and attachments that lie to 
the inside of the allowed attachments to for j ^ i. We include in this class attachments 
to the interior of eikonal lines. ("Interior" means to the inside of attachments of 
gluons.) 

The outermost allowed attachment of C* gluon to a singular line in {i ^ j) generally 
lies to the inside of attachments of additional gluons that have infinitesimally smaller energy 
scales. While the propagator immediately to the outside of the outermost allowed attach- 
ment of gluon is not precisely on the mass shell, it is on the mass shell, up to relatively 
infinitesimal corrections. Furthermore, if it is a gluon propagator, then its polarization is 
orthogonal to its momentum, up to relatively infinitesimal corrections. Therefore, when we 
apply the C"*" decoupling relation, no eikonal- line contribution appears at this point. 

The result of the application of the decoupling relations to the large-scale singular 
gluons with the largest energies is that the connections of these gluons to subdiagrams other 
than J* and the 6-quark line are replaced with connections to eikonal fines. The 
eikonal lines attach to the external-fermion lines just to the outside of H. (Here, and in 
subsequent discussions, "external-fermion lines" denote the fermion lines that originate in 
the external quarks and antiquarks that are associated with the mesons in our model.) The 
C' eikonal lines attach to the external 6-quark line and the external light-quark line from 
the B meson just to the outside of H. 

We can iterate this procedure for large-scale singular gluons with successively lower 
energy scales. After each iteration, there is a new C* eikonal line that attaches to each C* 
external-fermion line just to the inside of the C* eikonal line from the previous iteration. It 
is easy to see that, for each external-fermion line, the new eikonal line can be combined with 
the eikonal line from the previous iteration to form a single eikonal line, on which the 
singular gluons with lower energy scale attach to the outside of the C singular gluons with 
higher energy scales. Other orderings of the attachments yield vanishing contributions. We 
continue iteratively in this fashion until we have factored all of the large-scale C" gluons. 
After this decoupling step, the sum of diagrams represented by Fig. 4 becomes a sum of 
diagrams represented by Fig. 5. 
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FIG. 5: Diagram representing the sum of diagrams that occurs after one apphes the decouphng of 
the large-scale collinear gluons that is described in Sec. IIIH 1. 

2. Initial factorization of the nominal-scale C* gluons 

Next, we factor the nominal-scale C* singular gluons. In addition to the attachments 
enumerated in the case of the large-scale C* gluons, we include attachments to the nominal- 
scale S gluons. Then, the application of the decoupling relations leads to eikonal lines 
that attach to the following locations: to the external-fermion lines just to the inside 
of the large-scale eikonal lines from the previous step; to the nominal-scale S singular 
gluon lines just to the inside of the connections of those lines to the external-fermion 
lines. After this decoupling step, the sum of diagrams represented by Fig. 5 becomes a sum 
of diagrams represented by Fig. 6. Application of the C' decoupling relation leads to C' 
eikonal lines that attach to the following locations: to the external-fermion lines from the 
B meson just to the inside of the large-scale eikonal lines from the previous step; to the 
nominal-scale S singular gluon lines just to the inside of the connections of those lines to 
the external-fermion lines from the B meson. The C' eikonal line that attaches to a given 
external-fermion line from the B meson can be combined with the large-scale C' eikonal line 
from the previous step to form a single C' eikonal line. 
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FIG. 6: Diagram representing the sum of diagrams that occurs after one apphes the initial decou- 
pHng of the nominal-scale collinear gluons that is described in Sec. IIIH2. 

3. Factorization of the nominal-scale S gluons 

We now wish to apply the soft decoupUng relations to factor the nominal- scale soft gluons. 
In order to do this, we implement the approximations for the allowed attachments of 
the soft gluons to J^. (Recall that we do not apply the soft approximations or the soft 
decoupling relations to the attachments of the soft gluons to the external 6-quark line or to 
JK) On the connections to the subdiagrams, we modify the soft approximation in the 
following way: We combine the momentum of the nominal-scale soft gluon with the total 
momentum of the attached nominal-scale eikonal line from the previous step. Then, 
when we implement the 5*^ decoupling relations, the nominal-scale eikonal lines are 
carried along with the nominal-scale soft-gluon attachments. We apply the decoupling 
relations to the allowed attachments of the soft gluons to J^. We also include vanishing 
connections of the nominal-scale soft gluons to the interior of the large-scale eikonal 
lines [25]. The propagator that lies to the outside of the outermost allowed connection of a 
nominal-scale soft gluon to a line in is on shell, up to relative corrections of infinitesimal 
size. Furthermore, if it is a gluon propagator, its polarization is transverse to its momentum, 
up to relative corrections of infinitesimal size. Therefore, when we apply the 5*^ decoupling 
relations, no S"^ eikonal lines appear at those points. 
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FIG. 7: Diagram representing the sum of diagrams that occurs after one apphes the decouphng of 
the nominal-scale soft gluons that is described in Sec. Ill H 3. 

The result of applying the decoupling relations is that soft gluons attach to eikonal 
lines, to the external 6-quark line and to J'. The eikonal lines attach to the external- 
fermion lines just to the outside of the nominal- scale eikonal lines and just to the inside 
of the large-scale C"^ eikonal lines. Associated with each connection of a nominal- scale 
soft gluon to an eikonal line is a eikonal line. Associated with each connection of 
a nominal-scale soft gluon to the external 6-quark line or to is a eikonal line. Our 
sample diagram is now given by Fig. 7. 

4- Further factorization of the nominal- scale gluons 

We next factor the nominal-scale gluons from the S"^ eikonal lines. In order do this, 
we include formally the vanishing contributions that arise when one connects the nominal- 
scale gluons to all points on the S"^ eikonal lines that lie to the inside of the outermost 
connection of the nominal-scale soft gluons. We also make use of the following facts: a 
nominal-scale eikonal line that attaches to one of the external-fermion lines is identi- 
cal to the eikonal line that one would obtain by applying the decoupling relation to the 
attachments of the nominal-scale C"^ gluons to an on-shell fermion line (that does not have 
exactly momentum); a nominal-scale eikonal line that attaches to a nominal-scale 
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FIG. 8: Diagram representing the sum of diagrams that occurs after one apphes the further de- 
couphng of the nominal-scale collinear gluons that is described in Sec. IIIH4. 

gluon is identical to the eikonal line that one would obtain by applying the decoupling 
relation to the attachments of nominal- scale gluons to an on-shell nominal-scale soft- 
gluon line. Then, applying the decoupling relation, we find that the nominal-scale 
gluons attach to eikonal lines that attach to the external-fermion lines just to the inside 
of the large-scale eikonal lines. This situation is represented by the diagram that is 
shown in Fig. 8. 

The nominal-scale eikonal lines can then be combined with the large-scale eikonal 
lines. After performing those steps we arrive at the final factorized form for our sample 
diagram, which is given in Fig. 9. 

5. Completion of the factorization 

Now we can iterate the procedure that we have given in Sees. IIIH1-IIIH4, taking 
the nominal scale to be the next smaller soft-gluon scale. In these subsequent iterations, 
we include the connections of soft and collinear gluons that have already been described. 
In addition, we include formally, in the steps of Sees. Ill HI and IIIH2, the vanishing 
contributions from the connections of the large-scale and nominal-scale gluons to the 
soft gluons of higher energies and to the eikonal lines that are associated with those soft 
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FIG. 9: Diagram representing the sum of diagrams that occurs after one completely decouples the 
large-scale collinear gluons and the nominal-scale soft and collinear gluons. 

gluons. 

Proceeding iteratively through all of the soft-gluon scales, we produce new nominal-scale 
eikonal lines at each step that connect to the external fermion lines just to the outside 
of the existing S"^ eikonal lines. Each gluon that attaches to a nominal- scale S"^ eikonal line 
has attached to it a eikonal line. In addition, there are nominal-scale eikonal lines 
from the steps of Sec. Ill 112 that attach to the external-fermion lines just to the inside 
of the nominal-scale eikonal lines. After the further factorization of the nominal-scale 
gluons that is described in Sec. IIIH4, both of the S"^ eikonal lines that attach to a 
given external-fermion line can be combined into a single eikonal line. 

At each step in the iteration, new C' eikonal lines appear that attach to the external- 
fermion lines from the B meson just to the inside of the C' eikonal lines from the previous 
step. For each external-fermion line, the new C' eikonal line can be combined with the C' 
eikonal line from the previous step to form a single eikonal line. Similarly, at each step in 
the iteration, new C' eikonal lines appear that attach to the nominal-scale S singular gluon 
lines that attach to the external fermion lines from the B meson. These new C' eikonal lines 
attach just to the inside of the C' eikonal lines from the previous iteration. Again, for each 
external-fermion line, the new C' eikonal line can be combined with the previous C' eikonal 
line to form a single C' eikonal line. 
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FIG. 10: Illustration of the factorization for the case of e'^e" annihilation. After the use of the 
decoupling relations, gluons with momenta in the S singular region attach to eikonal lines and 
gluons with momenta in the singular regions attach to eikonal lines. 

Following this procedure, we arrive at the factorized form for the singular contributions. 
The S subdiagram now connects only to eikonal lines, to the external 6-quark line, and 
to JK The S*^ eikonal lines attach to the external-fermion lines just outside of H. All of 
the C"^ singular contributions are contained in the subdiagram and the associated 
eikonal lines, which attach to the external-fermion lines just outside of the 5*^ eikonal 
lines. All of the C' contributions are contained in the J' subdiagram and associated C' 
eikonal lines. These eikonal lines attach to the external-fermion lines from the B meson 
just to the outside of the H subdiagram and to S singular gluon lines just to the inside of 
the connections of those lines to the 6-quark line and to J'. This factorization is illustrated, 
for the case of e~^e~ annihilation, in Fig. 10. 

I. Forms of the S and functions and cancellations of eikonal lines 

1. Cancellations of the soft eikonal lines 

At this point, in the case of e'^e~ annihilation into two quarkonia, the 5* subdiagram and 
associated soft eikonal lines, which we call S, take the form of the vacuum-expectation value 
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of a time-ordered product of four eikonal lines: 

S{xig,Xiq,X2q,X2g) = (0|r{ [xi^, 00+] [oO+ , Xig] (g) [X2q, 00~][oO~ , X2q]}\0) S , (41) 

where Xiq and Xiq are the points at which the eikonal lines attach to the quark and antiquark 
external lines from meson 1 and meson 2, respectively, [y, x] is the eikonal line that is defined 
in Eq. (6), oo+ = (oo, 0, Oj_), and oo~ = (0, oo, Oj_). The symbol (8) indicates a direct product 
of the color factors that are associated with the soft-gluon attachments to meson 1 and the 
soft-gluon attachments to meson 2. The S subscript on the matrix element indicates that 
only the contributions from the S singular region arc kept. 

In the case of i?-meson decays, the S subdiagram is still connected to the B meson and 
takes the form 

SB{Xlq,Xiq,X2q,X2q) = (0|T{[xig, 00+] [oO+, Xig] (g) [X2q, 00~][oO~ , X2q\}'^ iT Bm^ b\B) ^ . (42) 

Here, we have suppressed the eikonal lines that are associated with the B meson. No 
soft gluons attach to those lines. 

Because the H subdiagram is insensitive to a momentum in the S singular region that 
flows through it, one can ignore the difference between xig and Xig and the difference between 
X2g and X2g- Therefore, in consequence of the fact that the external mesons are color singlets, 
the S'^ quark and antiquark eikonal lines cancel, and the S~ quark and antiquark eikonal 
lines cancel. In the case of e'^e~ annihilation into two quarkonia, this cancellation implies 
that the S subdiagram is completely disconnected, and, therefore, can be ignored. In the 
case of S-meson decays, the remaining S subdiagram now connects only to the external 
6-quark line and to J'. 

2. Rearrangement of the B-meson singular contributions 

As we have noted, there are C' eikonal lines associated with the B meson. These 
eikonal lines attach to the external-fermion lines from the B meson just to the outside of 
the H subdiagram and to S singular gluon lines just to the inside of the connections of 
those lines to the external 6-quark line and to J'. We can now remove the latter class of 
eikonal hues as follows. We note that, because the S subdiagram now connects only to the 
external 6-quark line and to J\ the C' eikonal lines that attach to the S singular gluons 
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are precisely the eikonal lines that would appear if one were to factor connections of 
singular gluons from S. (One can carry out the factorization iteratively, level-by-level, 
factoring the nominal-scale C' gluons from the nominal-scale S singular gluons.) Therefore, 
we restore the connections of the singular gluons to S and drop the C' eikonal lines that 
attach to S singular gluons. 

3. Forms of the meson distributions 

We make a Fierz rearrangement to decouple the color structures of the subdiagrams, 
the 6-quark and J' subdiagram, and their associated coUinear eikonal lines. Then, these 
subdiagrams and their eikonal lines are given by the following matrix elements: 

/+0O 
d\2x)exp[-2iq^ ■ x]{H,(P,)\^^(x)T{[x,oo-][oo-,-x]}^^(-x)\0)c+ (43) 
-oo 

for the C"*" quarkonium, 

/ + 00 
d\2x)eM-m2-x]{H2{P2)\^a{x)T{[x,oo+][oo+,-x]}^^{-x)\0)c- (44) 
-oo 

for the C~ quarkonium, 

/ + 00 
d\2x)exp[-2tfk-x]{K{pK)\^^{x)T{[x,^+][<^+,-x]}^^{-x)\0)c- (45) 
-oo 

for the light meson, and 

/+0O 
d^xeMm-xm'^i/s{x)T{[x,oo-][<x,-,0]}'^bamB{pB))s,ci- (46) 
■oo 

for the B meson. In Eqs. (43)-(46), a and /3 are Dirac indices. It is understood that the 
fields ^ and ^ in each matrix element are in a color-singlet state. In these distributions, 
the arguments qi, q2, and fk are each half the difference between the quark momentum 
and the antiquark momentum at the points at which they enter H, and the argument qi 
is the antiquark momentum at the point at which it enters H. We have suppressed the 
dependences on the total meson momenta Pi, P2, Pk, and in the arguments on the left 
sides of Eqs. (43)-(46). The subscripts S, , C~ and C' on the matrix elements indicate 
that we are retaining only the C'^, C~, and C' singular contributions. 
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4- Light-cone distributions and cancellations of the collinear eikonal lines 



Away from the endpoint region, we can simplify the factorized expression further. 

In H, away from the endpoint region, we can approximate the momenta of the quark and 
antiquark in the hght meson by their minus components. The leading relative errors in this 
approximation are of order qu/pK ~ Aqcd/Q- Then, integrating Jk over and ffc±, we 
obtain 



= ^ / dx+exp[-i(2y - 

TT J-oo 



x{K{pK)\^a{x+)T{[x+,oo+][^+,-x+]}^0{-x+)\O)c-. (47) 

The quantity H has been analyzed in the context of soft-coUinear effective theory (SCET) 
for the case of 5-meson decays into a lepton pair plus a photon [31] and for the contribution 
to 5- meson decays into two light mesons that arises away from the endpoint region [22]. 
The conclusion of these analyses is that H is given, to leading order in Aqcd/Q, by a matrix 
element of a SCET operator that depends only on the plus component of the momentum of 
the light antiquark in the B meson. ^ Furthermore, the SCET operator has a Dirac-matrix 
structure such that only the S-meson light-cone distribution contributes. We assume 
that a similar SCET analysis holds in the case of B meson decays to a quarkonium plus 
a light meson away from the endpoint region. Then, integrating Jb over pj' and pi±, we 
obtain 



J _ Pb dpi d'Pi^ - 



2n 



-oo 



dx exp[i^p^a; ] 
x{0\^if,{x-)T{[x-,oo-][oo-,0]}^bamB{pB))c+, (48) 



where ^ = p^ /Q. 

We do not approximate the momenta of the heavy-quark and heavy antiquark in the 
quarkonia by their dominant momentum components because, in so doing, we would intro- 



^ These analyses arc based on Lorentz (rcparametrization) invariance and power counting in -^/Aqcd/Q- 
The next-to-leading-order spectator-scattering contributions to S-meson decays to light mesons have been 
computed in Refs. [32-36] and confirm the general analysis for this process. 
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duce errors of relative order rricv/Q for each quarkonium. As we will explain in Sec. IIIK, 
such an error would be larger than the errors that arise from the approximations that we 
have used to derive the factorization result. 

Now, we can see that there is a partial cancellation of the C~ quark and antiquark 
eikonal lines in Eq. (47) and a partial cancellation of the quark and antiquark eikonal 
lines in Eq. (48). The cancellations would be complete, were it not for the fact that the 
H subdiagram is sensitive the routing of coUinear momenta through it. This sensitivity 
corresponds to the separation in space-time of the points x'^ and —x'^ in Eq. (47) and the 
points x~ and in Eq. (48). The quark and antiquark eikonal hues in Eqs. (47) and (48) 
cancel where they overlap, leaving an eikonal line that runs directly between the quark and 
the antiquark: 

JKap{y) = ^ r"cia;+expH(2y-l)p^a;+](ir(p^)|^r,(a;+)P[a;+,-a;+]^^(-a;+)|0)c- 

J-oo 

= Y.^KAyWKj],p. (49a) 

3 

JfiaMO^^ r rfx-exp[iepSx-](0|*,^(x-)P[x-,0]*;.„(0)|S(pB))c+ 

= Y.^Bm{0[TBmU. (49b) 

m 

where we have written the time-ordered product of the exponentiated line integral as a 
path-ordered product.^ The expressions in Eqs. (49) have the form of the conventional 
light-meson and S-meson light-cone distributions, but, at this stage, they contain only the 
singular contributions to those light-cone distributions. Since the integrations over y and 
^ have a finite range of support in the typical separation of the points x+ and — x+ in 
Eq. (49a) and the points x~ and in Eq. (49b) is of order 1/Q. 



J. Factorized form 



1. Factorization of the logarithmic enhancements 

At this point, we have established that the contributions from the soft singular region 
decouple completely from the subdiagrams (leaving no residual eikonal lines). We have 

^ Reference [15] contains an incorrect statement that the eikonal lines in Eq. (49) cancel completely. 
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also established that the contributions from the coUinear singular regions factor from the 
H subdiagram and are contained entirely in the J^, Jk and Jb subdiagrams. As we have 
mentioned, in the case of e'^e" annihilation into two quarkonia, the S subdiagram is now 
completely disconnected, and can be ignored. In the case of 5-meson decays, the S subdi- 
agram is still connected to the external 6-quark line and to J/ {i.e., to Jb). 

Now let us restore rric to its nonzero physical value. Then, some of the soft and coUinear 
singularities become would-be soft and coUinear singularities. However, the would-be sin- 
gularities are still contained in the J^, Jk and Jb subdiagrams. Therefore, there are no 
actual or would-be coUinear singularities in the H subdiagram. Furthermore, there are no 
actual or would-be soft singularities in the H subdiagram. In the case of 5-meson decays, 
there are, however, soft contributions from the endpoint region in the S — S subdiagram, 
and, hence, in the H subdiagram. As we have emphasized, these endpoint contributions are 
associated with the topology of Fig. 2(b). 

Next let us redefine J^, Jk and Jb by extending the ranges of integration from the 
infinitesimal C^, C' singular regions and, in the case of Jb, the S singular region, to finite 
regions that are defined by an ultraviolet cutoff /ii? ~ Q on the logarithmic integrals. H is 
then redefined to be the remainder of the amplitude. One can think of fip as an infrared 
cutoff on the soft and coUinear enhancements in H. This redefinition has the effect of 
absorbing the coUinear logarithmic enhancements that are associated with the coUinear 
singularities into J^, Jk and Jb- It also has the effect of absorbing soft enhancements that 
are associated with soft singularities into Jb- 

One might worry that, in making such an extension, we could introduce new singularities 
and logarithmic enhancements in J^, Jk and Jb that are associated with their coUinear 
eikonal lines. The lightlike eikonal lines that are parametrized by the vectors ni, n2, and 
fii could, in principle, be sources of gluons that are coUinear to the minus, plus, and e/ 
directions, respectively, as well as sources of soft gluons. In fact, this does not happen in 
the case of the hght-meson light-cone distribution [Eq. (49a)] or the S-meson hght-cone 
distribution [Eq. (49b)]. As we have noted, there is a partial cancellation between the quark 
and antiquark eikonal lines in these light-cone distributions. The remaining cikonal-line 
segment is typically of length 1/Q. Therefore, only modes with virtuality of order Q can 
propagate along it, and no coUinear or soft singularities or logarithmic enhancements are 
associated with it. 
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In the case of the distributions in Eqs. (43) and (44) and the Jk and Jb distributions in 
Eqs. (45) and (46), which are appropriate when the hght-meson momentum is in the endpoint 
region, we make use of a trick to prevent coUinear singularities and enhancements from 
developing along the eikonal lines: In each case, we replace the lightlike eikonal lines with 
spacelike eikonal lines. That is, we replace the eikonal-line vectors ni, and fii with a vector 
Uz — (l/-\/2)(l, —1, 0_l), which points in the z direction. Because of the freedom in choosing 
the coUinear eikonal vectors that we described in Sec. Ill E, this replacement has no effect on 
the C"^, C~, and C' singular contributions in J+, J", Jk, and Jb, respectively. Furthermore, 
the soft singularities (and enhancements) that arise from soft-gluon attachments to the 
quark and antiquark eikonal lines in Eqs. (43), (44), (45), and (46) cancel. This cancellation 
derives from the following facts: The ^'-singular attachments lie to the exterior of any non- 
S'-singular attachments to the eikonal lines; any non-S'-singular attachments are within 1/Q 
of the eikonal-line endpoints; the endpoints —x and x in Eqs. (43), (44), and (45) and 
and X in Eq. (46) are within l/Q of each other. Hence, one can argue, as in Sec. Ill II, 
that the segments of the quark and antiquark eikonal lines that contain 5'-singular-gluon 
attachments cancel. 

We have argued that there arc neither soft nor coUinear logarithmic enhancements in 
the H subdiagram. Therefore, in the cases of e^e~ annihilation and i?-meson decay in the 
topology of Fig. 2(a), the H subdiagram involves only momenta of order Q. The lower- 
virtuality momenta are contained in the distributions in Eqs. (43) and (44), in 
Eq. (49a), and $b in Eq. (49b). 

2. Further factorization of the endpoint contributions 

In the case of i?-meson decays in the topology of Fig. 2(b), the H subdiagram is also free 
of soft and coUinear logarithmic enhancements, but it still contains gluons with momenta 
of order Aqcd that arise from the endpoint region. These gluons consist of the gluon that 
is marked with an asterisk in Fig. 2(b) and gluons that are radiated from it. They are the 
part ol S — S that remains after S has been extended to include soft enhancements. They 
can connect to active-quark or active- antiquark lines (those that participate in the weak 
interaction). However, they cannot connect to any part of the or J~ subdiagrams, which 
reside to the outside of the connections of the soft gluons to the active-quark or active- 
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antiquark lines. Because these soft gluons connect the B meson and hght meson to the 
quarkonia, they potentially violate the factorized form in the second term of Eq. (2). 

However, we can make a further decoupling of the connections of the endpoint soft gluons 
from the active-quark and active-antiquark lines in the quarkonium. We apply a modified 
soft approximation to these gluons. Because the soft gluons have a finite soft momentum 
of order Aqcd, rather than a soft singular momentum, there are errors associated with 
the application of the soft approximation to the quark or antiquark line that are order 
Aqcd/Q- These errors are negligible in comparison with the errors that are associated 
with the modified soft approximation for the average of the quark and antiquark momenta. 
Next we apply the decoupling relation. Then, the soft gluons attach to eikonal lines 
that attach to the heavy-quark and heavy- antiquark lines just outside the H subdiagram. 
Because the remaining part of H is insensitive to routing of the soft momenta through it, 
the quark and antiquark eikonal lines cancel, up to corrections of order Aqcd/Q- Then, the 
endpoint contributions are contained entirely in a subdiagram BK, which consists oi S — S 
(after S has been extended to include soft enhancements), the parts of the S- meson and 
light-meson quark and antiquark lines to which S — S attaches, Jk in Eq- (45), and Jb in 
Eq. (46).. The H subdiagram now contains only momenta of order Q. Consequently, we 
can contract H to a point with respect to the soft interactions in BK. Then, decoupling 
the Dirac and color indices of BK from H by making Fierz rearrangements, we obtain the 
S-meson-to-light-meson form factors in Eq. (3) from BK and short-distance coefficients He 
from H. 

3. NRQCD decomposition of the quarkonium distribution amplitudes 

At this stage, we have achieved the factorized forms of Eqs. (1) and (2), except that the 
quarkonium factors are expressed in terms of quarkonium distribution amplitudes, instead 
of NRQCD matrix elements. We now argue that the quarkonium distribution amplitudes 
can be expanded as a sum of products of NRQCD matrix elements times short-distance 
coefficients. 

The distribution amplitude describes the local creation of a quark- antiquark pair, 
followed by its evolution, through QCD interactions, into a quarkonium. The gluons in 
the distribution amplitude, which have momentum in the e'^e~ CM frame or the 



46 



S-meson rest frame, have hard, soft, and threshold (potential) momenta in the quarkonium 



and any number of gluons and light-quark-antiquark pairs, then it is clear that it could be 
written as a standard NRQCD decomposition of a full QCD amplitude. That is, it could 
be written as a sum over products of short-distance coefficients times matrix elements of 
local NRQCD operators. Lines with virtualities of order rric or greater lie to the inside of the 
lower virtuality lines, and could be integrated out to yield the local NRQCD operators times 
short-distance coefficients. Lines with virtualities less than of order rric are well described 
by NRQCD and would be accounted for by the NRQCD matrix elements of these local 
operators between the vacuum state and the quarkonium state. 

A complication to this picture arises because the distribution amplitudes also contain 
eikonal lines, which are not a part of QCD. However, the attachments of the eikonal lines to 
the external heavy-quark and heavy- antiquark lines at the points Xig and Xig are separated 
in space-time by a distance of order 1/Q. Hence, only high- virtuality modes can propagate 
on these hues. Therefore, they too can be integrated out to yield local operators times 
short-distance coefficients. These operators would involve the gauge field, as well as the 
quark and antiquark fields. 

Therefore, we can write 



where an and are short-distance coefficients, an and a2i each have two Dirac indices, 
corresponding to the quark line and the antiquark line in J+ and J~ , respectively. We 
suppress those indices. We then make the following identifications for the cases of e+e~ 
annihilation, B decay in the topology of Fig. 2(a), and B decay in the topology of Fig. 2(b), 
respectively: 



rest frame. If the distribution amplitude involved only a heavy quark, a heavy antiquark. 




(50a) 




(50b) 




(51b) 



(51a) 



(51c) 
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Here, we have also suppressed the Dirac indices on H, which are contracted into the (sup- 
pressed) Dirac indices on an and a2j- The identifications in Eq. (51) lead directly to the 
factorization formulas in Eqs. (1) and (2). 

K. Corrections to factorization 

Now let us discuss the corrections to the factorized form. The most important corrections 
to the factorized form arise because of the approximate nature of the cancellations of the 
couplings of soft-singular gluons to the color-singlet quarkonia. These cancellations hold 
only up to the errors in the modified soft approximation. We wish to compare the sizes of 
these errors relative to the factorized contributions. In some cases, the contributions from 
the modified soft approximation, which ultimately cancel, simply scale with mc, Q, and v 
in the same way as the factorized contribution. However, there can be exceptions to this 
scaling because of the specific quantum numbers of the final states in a given process. We 
give some examples of such exceptions below. 

Note that, in the case of e'^e" annihilation into two quarkonia, violations of factorization 
arise only in contributions involving the corrections to the modified soft approximation for 
both quarkonia. The reason for this is that, if the S subdiagram decouples from quarkonium i, 
but not from quarkonium j, then the S subdiagram can be absorbed into the definition of 
the J subdiagram for meson j. (See Ref. [16] for a more detailed discussion of this point.) 

Now let us discuss the dependence of the relative size of the corrections to factorization on 
the orbital angular momenta of the produced quarkonia. The leading errors in the modified 
soft approximation are proportional to qi±/Pi ~ rricv/Q. Because of their proportionality 
to qj_L, the leading errors in the modified soft approximation contribute one unit of orbital 
angular momentum. Consequently, in order to yield a QQ pair in the quarkonium angular- 
momentum state, they must be accompanied by an additional factor qi±/'mc ~ v from the 
short-distance production process in the case of an 5'-wave quarkonium and an additional 
factor {qi±/mc)^~^ ~ v'"^^ from the short-distance production process in the case of an 
L-wave quarkonium with L > 0. The factorized contributions contain a factor for each 
L-wave quarkonium. Therefore, the factorization-violating contributions arc suppressed, 
relative to the factorized contributions, by a factor fi ~ rricv'^/Q for each quarkonium i in 
an 5"- wave state and by a factor /j ~ nic/ Q for each quarkonium i in a higher orbital-angular- 



48 



momentum state.^ The suppressions of the factorization- violating contributions that we find 
here are consistent with those that were found in Ref. [15]. However, in Ref. [15], powers of 
V in the factorization-violating contributions were ignored. 

The relative sizes of the corrections to factorization can depend on additional quantum 
numbers, beyond the orbital angular momenta of the quarkonia. Let us mention a few ex- 
amples. In the case of production of S'-wave quarkonia, the factorized production process 
can be suppressed by powers of rric/Q if it involves a helicity flip. (See, for example Ref. [8].) 
However, we expect such a helicity suppression to apply to the factorization-violating con- 
tributions, as well, and so it should not affect the relative size of the factorization- violating 
contributions. In order i?-meson decays do not produce a Xco oi Xc2 charmonium (a 
J = or J = 2 P-wavc state). Those processes are allowed only in order a^. On the 
other hand, the factorization-violating corrections to S-meson decays do produce XcO and 
Xc2 charmonia in order q;°. Therefore, the factorized process for Xcq or Xci production is 
suppressed by a power of cts, relative to the factorization- violating process, and may not be 
dominant. Since the factorization- violating contributions arise from diagrams in which at 
least one gluon has been added to the leading-order process, there can also be a dependence 
of the relative size of the factorization- violating contributions on the color structure of the 
hard subprocess. 

In perturbation theory, the factorization-violating contributions may be enhanced by 
logarithms of /ml. Furthermore, they are infrared divergent. In reality, these infrared 
divergences are cut off by nonperturbative effects associated with confinement. Our analysis 
does not determine the size of these factorization- violating contributions: It only shows that 
they vanish as one or two powers of / as / approaches zero. One might use the small 
parameter / as an estimate of the size of the factorization- violating contributions. However, 
the size of the factorization- violating contributions is an issue that, at present, must be 

^ In the case of the factorized form for _B-meson decays in the first term of Eq. (2), there are also errors 
in the cancellation of the couplings of the soft-singular gluons to the light meson. These errors are of 
order qk/Q ^ Aqcd/Q relative to the factorized contributions. In addition, there are errors of relative 
order Aqcd / Q that arise when one expresses the amplitude in terms of the light-cone distributions for the 
light meson [Eq. (47)] and B meson [Eq. (48)]. These errors arise because one neglects in H the plus and 
transverse components of the momenta of the quark and the antiquark in the light meson and the minus 
and transverse components of the momentum of the antiquark in the B meson. We neglect these errors 
in comparison with the errors in the cancellation of the couplings of soft-singular gluons to the quarkonia. 
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settled through experiment or, perhaps, lattice simulations. 

We note that, at lowest order in v, one sets qi — 0, and the cancellation of the couphngs 
of soft-singular gluons to each quarkonium is exact. At lowest order in v, only 5"- wave 
quarkonium production is possible. An explicit calculation of the one-loop corrections to 
S'-wave quarkonium production in S-meson decays at lowest order in v [13] confirms our 
expectation that these corrections are free of infrared divergences. In the case of double 
quarkonium production in e'^e^ annihilation, the exact cancellation of the couplings of soft- 
singular gluons holds for each quarkonium that is treated at lowest order in v. If only one 
quarkonium is treated at lowest order in v, then S can be absorbed into a re-definition of 
the distribution function of the remaining quarkonium, and the cancellation of factorization- 
violating infrared divergences is expected to be exact. This expectation is confirmed by an 
explicit calculation of the one-loop corrections to a[e'^e~ J/ip + Xcj], where the J/ip is 
treated at lowest order in v [37]. Even in the case of S-wave quarkonium production, we 
expect infrared divergence to appear at higher orders in v, accompanied by a suppression 
factor /, as discussed above. 

Finally, we mention that we could have written the coUinear functions that are asso- 
ciated with each quarkonium in terms of light-cone distributions instead of NRQCD matrix 
elements. The derivation of this result would entail the use of coUinear approximations for 
the momenta of the heavy-quark and heavy- antiquark in meson 1 (2) in which one neglects 
the minus (plus) and transverse components in comparison with the plus (minus) compo- 
nents. These approximations introduce an error of relative order fi for each quarkonium, 
and, in the case of double-quarkonium production, these errors must be added, rather than 
multiplied, in order to obtain the error for the complete amplitude. In the resulting fac- 
torized expression, no power-suppressed soft divergences would appear in H because, when 
the quark and antiquark momenta in each quarkonium are taken to be coUinear to each, 
the cancellation of soft divergences between the quark and antiquark in each quarkonium is 
exact. In contrast, in the factorized expressions involving NRQCD matrix elements that we 
have presented, power-suppressed soft divergences do appear in H and must be discarded 
in order to obtain the factorized expression. However, as we have said, these divergences 

This result falsifies the conjecture in Ref. [15] that the soft cancellation might be inexact in higher orders 
in Ug, even at lowest order in v. 
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are suppressed as fifj (rather than /j + fj) in double-quarkonium production. Hence, in 
the case of double-quarkonium production, a factorized expression involving light-cone dis- 
tributions for the quarkonia would be less accurate than the expression involving NRQCD 
matrix elements that we have presented. 



IV. ONE-LOOP EXAMPLES 



In this section, we illustrate some of the features of the factorization result by presenting 
some one-loop examples for double-charmonium production in e'^e" annihilation and for 
production of a charmonium and a light meson in S-meson decay. In our examples, we wish 
only to identify the soft divergences, and so we consider the loop gluon to be soft, and we 
make use of the soft approximation in our calculations. 



A. Soft approximation 

In the general factorization argument, we have taken the soft approximation to be the 
replacement 

in the gluon-propagator numerator. Recall that k is the soft-gluon momentum, that p is 
the momentum of the line in the coUinear subdiagram, and that /i and u correspond to 
the attachments of the gluon to the coUinear and soft subdiagrams, respectively. Consider 
now, for instance, an initial-state quark with momentum p that absorbs a gluon with soft 
momentum k. Then, we can make use of the graphical Ward identity (Feynman identity) 
to rewrite the soft approximation: 

k^tPu \ {:^ + ^) +m 



k-p + iej (p + kf - m2 + ie^^^^^'' 



k-p + ie(p + ky -m'^ + ie \k-p + ie^ 

(53) 

Since we are interested only in identifying the soft divergences, we eliminate any ultraviolet 
divergences in loop integrals, without affecting the soft divergences, by reintroducing the k^ 
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terms in the quark and antiquark denominators. That is, we make the substitution 

2k-p + ie ^ k"^ + 2k-p + ie (54) 

in the denominator of the last line of Eq. (53). This is the form of the soft approximation 
that we will use. Therefore, if a quark or antiquark line with physical momentum pi absorbs 
a gluon with soft momentum k, color a, and vector index /i, then the full amplitude is 
approximated as 

A[Q{p.) + gm ^ 9sT^ ^,^^J:^^^^^ A[Q{p.)], (55) 



(56a) 



(56b) 



where 

1 initial-state particle 
-1 final-state particle 

1 quark 
-1 antiquark 

Some of the calculations that we present below involve soft-gluon loop corrections in 
which a gluon can be emitted or absorbed by two particles, each of which can be a quark or 
an antiquark. [See Figs. 11 and 12(a)-(d)]. If we choose the sense of the loop momentum k 
such that it is absorbed by the line with momentum p^, then application of Eq. (55) yields 
the soft loop factor 

w \^_-2f_^ "^bjbjPi ■ Pj 

[Pi, Pj) ^9s J ^2j:)d [A;2 + 2aik ■ Pi + ie] [A;^ - 2ajk ■ pj + is] [A;^ + is] ' ^ ' 

where we have regulated the soft divergence by using dimensional regularization, with d — 
4 — 2e. The infrared-divergent part of the multiplicative correction factor from this soft loop 
factor is given by 

where pi is the physical momentum of the particle i, 



(See also Ref. [5]). 



2 2 

A,^^(p.,P,) = ^/l-^^^, (59a) 
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FIG. 11: One-gluon corrections to the double-quarkonium production amplitude. The blob labeled 
H represents the lowest-order hard-scattering process in which two heavy quark-antiquark pairs 
are created. 

B. Exclusive double quarkonium production 

In this section we consider double quarkonium production in e~^e~ annihilation, 7* — )■ 
Hi{Pi) + (-P2)- Hi and H2 are the produced quarkonium states. For definiteness, we will 
assume that the produced quarkonia are charmonium states. 

Recall that there are QCD and QED contributions to the Born-level hard-scattering pro- 
cess for double-quarkonium production in e^e^ annihilation. (See Ref. [8] for details.) How- 
ever, since this is an exclusive process, only color-singlet QQ pairs can contribute. Therefore, 
the soft-gluon loop corrections to the QED diagrams are zero at one loop. Nonvanishing 
corrections to the QED diagrams appear only at two-loop order. 

We now write the amplitude as 



where C is the appropriate color factor and Rhi+H2 is the soft loop factor, which is given by 



We obtain I{pii,p2j) from Eq. (58). Retaining only those terms with one or fewer factors of 



^soft 



(60) 



(61) 
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qi or q2 and discarding terms containing two or more powers oi z — Arn^/s, we obtain 



a. 



47reiR 



X 



2bibj'Ki + bibj (2 In z + Az) - 46, 



Pi-q2 
Pi-P2 



Pl-P2 



[Pi ■ P2 



(62) 



where 6j is defined in Eq. (56b). 

Now let us evaluate the various scalar products that appear in Eq. (62). Prom Eqs. (20), 
(22), and (23), we have 

1 



Pl-P2 

Pi-q2 
P2-qi 
qi ■ q2 



Q.lQ.2 



= -Pi ■ P2 



QtQ2 
M1M2 



M1M2 

qi± ■ q2±, 



+ qi± ■ q2± 



where 



CM 



(63a) 
(63b) 
(63c) 

(63d) 
(64) 



Thus, wc sec that all of the invariants are of order Q"^ = s and that the various terms in 
Eq. (62) are of order (mc/Q)^. If we add the contributions of Figs. 11(a) and 11(b) or 
Figs. 11(c) and 11(d), then the first, second, and third terms in Eq. (62) cancel. Similarly, 
if we add the contributions of Figs. 11(a) and 11(c) or Figs. 11(b) and 11(d), then the first, 
second, and fourth terms in Eq. (62) cancel. Therefore, we obtain 



Now, 



_ gi • g2 gi • -P2g2 • Pi 

""^^^^^ " TreiR [ Pi • P2 (Pi • P2)^ , 

2qi± • q2i. 



Pl-P2 

qi ■ P2q2 ■ Pi 
(Pi ■ P2V 



+ 



M1M2 s - Mf - Mi 



2 ' 



Therefore, we have 



R 



H1+H2 — 



MiM2[s - Ml - M^f 
qi± q2± . .... 2M1M2 



(65) 

(66a) 
(66b) 



vreiR IS- Mf 
16as qi± ■ q2± 



Ml 



+ qiq2- 



Ml 



M2 



-z ~ 



{rricvf 



(67) 
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where we have used the fact that the components of qi±_ are of order rricV. Equation (67) 
shows exphcitly the suppression of the soft divergence that is expected from the factorization 
proof. 

We can see that the soft divergent terms are proportional to one power of qi and one 
power of ^2- Therefore, double- S'- wave production (7* — )■ J/tp + r]c) and ^-wave/P-wave 
production (7* J/%1) -\- Xc, Vc + he) at one- loop and at leading order in v are free of soft 
divergences. This is confirmed by the exphcit one-loop calculations of Refs. [37, 38], for the 
J/ip + r]c and J/ip + Xc cases. 

Note that, if we consider only the contributions from the diagrams of Figs. 11(a) and 
11(b) or Figs. 11(c) and 11(d), then the fourth term in Eq. (62) survives. Similarly, if we 
consider only the contributions from the diagrams of Figs. 11(a) and 11(c) or Figs. 11(b) 
and 11(d), then the third term in Eq. (62) survives. Each of these terms are of order 
(rric/Q)^. At first sight, this seems puzzling, since we expect soft divergences to cancel 
up to terms of order rric/Q when we add the contributions from the connections of a soft 
gluon to the quark and antiquark in a quarkonium. The failure of that cancellation in 
the present case can be understood because the soft approximation that we have taken 
contains enhancements that arise when the gluon momentum is nearly collincar to cither of 
the quarkonia momenta. Actual coUinear divergences (logarithms of rric) are absent in the 
third and fourth terms in Eq. (62) because, in the one-loop case, they appear with equal 
strength in the contributions in which the gluon connects to the quark or the antiquark 
in a quarkonium. However, a residual finite piece of the coUinear enhancement survives 
because the soft approximations for the quark and antiquark lines are not equal when the 
gluon momentum is in the coUinear region. This failure of the soft cancellation when the soft 
function contains coUinear enhancements was noted in Ref. [16]. In the proof of factorization 
that we have given, such coUinear enhancements are removed from the soft function S and 
reside in the functions and associated eikonal lines. In our one-loop example, we have 
neglected the dependence of the hard function on the momentum of the gluon. Therefore, 
the C"*" eikonal lines that arise from the C"*" enhancements in the diagrams of Figs. 11(a) 
and 11(c) or Figs. 11(b) and 11(d) cancel, and the C~ eikonal lines that arise from the C~ 
enhancements in the diagrams of Figs. 11(a) and 11(b) or Figs. 11(c) and 11(d) cancel. The 
sum of all four diagrams in Fig. 11 is therefore free of coUinear enhancements that could 
spoil the soft cancellation and is in accord with the results of the factorization proof. 
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FIG. 12: One-gluon corrections to the lowest-order B decay amplitude. The black square represents 
the electroweak interaction. 

C. S-meson decays 

In this section we consider decay of a. B meson into a light meson plus a charmonium 
state. For definiteness, we take the light meson to be a meson. Therefore we have the 
process B{pb) — )• H{Pi) + K{j)k)- H is the produced charmonium state, which we take 
to be a ^Pj state. The one-soft-gluon corrections to the lowest-order decay amplitude are 
represented in Fig. 12. We refer to diagrams (a)-(d) as vertex corrections and to diagrams 
(e) and (f) as spectator contributions. 

1. Vertex corrections 

In discussing soft contributions to the vertex corrections, we assume that the light quark 
has a small mass m^, and we work in the limit — 0. In order to make contact with 
the results in Ref. [39], we neglect qu = —ypx + Pk, in comparison with px, and we neglect 
Qb ~ Aqcd in comparison with m^. We keep terms containing zero or one power of qi. 

For the vertex corrections in which the soft gluon attaches to the light-quark line 
[Figs. 12(c) and 12(d)], we obtain 



^K.PM.-,J-l-8 ^,,^^/_;^,,. ^1^8^^J. (68) 



56 



The upper sign corresponds to the quark with momentum pig, and the lower sign corresponds 
to the antiquark with momentum pig. The soft loop factor is then given by 



47reiR 



In 



2 2 



{ml 



Pb ■ qi 



+ 2m 



(69) 



4.mlf ) ' {ml -Ami) 

This expression contains coUinear divergences, which manifest themselves when we take 
rUs — 0. The divergences arise because the soft expressions contain contributions from C~ 
momentum. However, as we have explained in Sec. IV B, at one-loop order, the coUinear 
divergences cancel when one sums over the connections to the quark and the antiquark in 
the charmonium. Computing that sum, we obtain 



I {Pk^Plq) +1 {Pk^Plq) 



16a 



s Pb ■ qi 



16a. 



Pb ■ qi 



(70) 



47reiR {ml - Amf) 47reiR ml{l- z)' 

where, again, z = iml/mf. 

For the vertex corrections in which the soft gluon attaches to the 6-quark line [Figs. 12(a) 
and 12(b)], we obtain 

p{pb,Pig{q)^ = A(l±5i), 

with 

64:mlmlpb ■ qi 



^ ml - Ami _ 

ml + 4ml^ {ml + 4ml){ml — 4ml)'^ 
The corresponding soft loop factors are given by 



(71) 
(72) 



l{pb,Plq(g)) = 



etc 



47reiR 



Pi 



In 



+ 



Pi 
2Pi 



(73) 



l + PiJ 1 

Summing over both diagrams in Figs. 12(a) and 12(b), we obtain 



HPb,Plq)+I{Pb,Plq) = 



2as 5i 
47reiR /3i 



In 



/3i 



/3i 



2/3i 



16as Pb ■ qi 2z 



In 2; + 



P! 

l-z^ 



(74) 



47reiR ml {\ - zf V"" ' ' 2z 
As we have explained in Sec. IV B, it is necessary to add the contributions of all four 
diagrams in order obtain the suppression of the infrared-divergent terms because the soft ex- 
pressions contain C+ contributions that spoil the soft cancellation. Those (7+ contributions, 
which correspond to eikonal lines in the general factorization proof, cancel at one-loop 
order when one sums over the connections to the h quark and the light quark. When we add 
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the contributions of all four diagrams, i.e., Eqs. (70) and (74), the leading term does indeed 
cancel, and we obtain 



The remaining infrared-divergent terms are suppressed at least as ^ In ^. We generally expect 
a suppression of the factorization- violating contributions by only a factor y/z. However, as we 
have mentioned in Sec. Ill K, the suppression factor can become z for production of P-wave 
quarkonia if one neglects the transverse momenta of the constituents of the B meson and the 
light meson, as wc are doing in the present example. As was noted in Ref. [39], the expression 
in Eq. (75) gives a nonvanishing contribution to the production of a ^Pj charmonium only 
if J = or J = 2. The Born- level cross section to produce a ^Pj charmonium with J = or 
J — 2 vanishes, and so the violations of factorization are suppressed only as z In z/ as with 
respect to the leading factorizing terms. 

Finally, we mention that, when we include the Born factors in the amplitude, along with 
the soft factor, and decompose qi and the quarkonium spin polarization e* into the J = and 
J = 2 angular-momentum tensors, then we obtain agreement with the results in Eqs. (14) 
and (16) of Ref . [39]. 

2. Spectator contributions 

In the spectator contributions of Figs. 12(e) and 12(f), we initially assume that pi, 
and ypK are all of order Aqcd- Then, because the gluon in Figs. 12(e) and 12(f) carries 
momentum p^- — pi ^ Aqcdj the heavy-quark or heavy-antiquark propagator is off shell by 
order mfeAqcD- That is, its momentum is in the semihard region. Therefore, in this example, 
we are discussing the further factorization of the gluons with endpoint soft momenta that was 
described in Sec. IIIJ2. Eventually, we wish to make contact with the results in Ref. [39]. 
In that work, q^ — ypx —Pkg was neglected in comparison with ypx- Neglecting generates 
endpoint divergences in y that are cut off in our model when ypx ~ Q'fc ~ Aqcd- Our 
discussion in this example also applies to the cancellation of those endpoint divergences. 

Since the gluon in Figs. 12(e) and 12(f) carries momentum of order Aqcd or less, we 
can apply the soft approximation to the connections of the gluon to the heavy-quark and 
heavy-antiquark lines. Keeping terms up to order qi, we find that the soft factor for the 



I {Pk„Plq) + I {Pk„Plq) + I {Pb,Plq) + I iPb,Plq) 




(75) 
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heavy quark and antiquark lines for sum of the diagrams in Figs. 12(e) and 12(f) is 

1 4 

where 

g=gi-Pi— -, 77 

Pi ■ [Pk, - Pi) 

and we have included the factor l/{pkq — Pi)^ from the gluon propagator in S*^. Wc note 
that Sfj_ is proportional to g ~ rricV. We will concern ourselves only with the contributions 
to the production a P-wave quarkonium at leading order in v. Therefore, we can retain 
only terms of leading order in v in the remaining factors that are associated with the heavy 
quarkonium. 

In computing the remaining factors in the spectator amplitudes, we make use of the 
quark-antiquark spin-projection operators that are given in the Appendix. In the projector 
for heavy quark and antiquark, we take the spin-triplet case, which corresponds to the 
calculation for the XcJ in Ref. [39]. The projector for a charm-quark pair production can 



be obtained by setting niq — rriq — rric and Eq + Eq — 1yjm\ -|- q| in Eq. (A9b). Retaining 
only the terms of leading order in v and using relativistic normalization, we have 

n~^-^/*(fi + 2me). (78) 

We obtain the P-meson projector by setting — m^, Eq — rrih, rriq — mi, and Eq — Ei in 
Eq. (A6) and retaining the terms of leading order in Aqcd- Using relativistic normalization, 
we have 

nf«CB(A + ^6)75(1-^), (79) 

where the light-antiquark mass mi is of order Aqcd and the factor Cb is defined by 

1 I mi 



^'-2\lm,{l + Ei/miy 
Similarly, we obtain the X-meson projector by retaining the terms in Eq. (A9a) of leading 
order in Aqcd- Using relativistic normalization, we have 



nf fti Ck 



mq 



^Kl5 - —i>Kiklb \ , (81) 



rriq 



y 

where we have retained small masses mq and m^ of order Aqcd for the quark and antiquark, 
respectively. The coefficient Ck is defined by 

mq{Eq + mq + Eg + mq) 
2^2{Eq + mq){Eq + mq){Eq + Eq) ■ 
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Now, the trace over the heavy quark and antiquark hnes is 

U, = Tr[n~7p(l - 75)] ^ -2x/2mee*. (83) 

The factor 7^(1 — 75) comes from the V — A weak vertex. The trace over the S-meson and 
hght-meson quark and antiquark hnes is 

L'; = Tr [nf 7. (1 - 75) nf 7''] 



CbCkTt {^Klb - ■^i>Kiklb^ Ip (1 - 75) (A + rn^) 75 - 



(84) 



where the factor 7^(1 — 75) comes from the V — A weak vertex, and, for simphcity, we have set 
'f^q — f^q — mi — m. The complete amphtude corresponding to the diagrams in Figs. 12(e) 
and 12(f) is 



^spect; 



ator 



3/2 



-s.u'^L'; 



-ig'^Cj.CnC'KCEw 



8\/2///, 



xTr 



-^Kj5 + 



{Pl - PkgfPl ■ {Pkg - Pi) 



K 



m 



/*(l-75)(A + m6)75 1-- 



m 



(85) 



Here, Cew = (G^f/"\/2)[K&Ks^i ~ ^6^s(C4 + Cq)], where Gf is the Fermi constant, the 
Vq^q^ are the CKM matrix elements, and the are the Wilson coefficients of the effective 
electroweak Hamiltonian. (See, for example, Ref. [39] for details.) 

In Eq. (85), the terms proportional to the spatial components of pi vanish upon integration 
over the angles that are associated with those spatial components. Then, using the fact that 
7o ~ terms of relative order Aqcb/ttii,, we can write 

-ig-CyC'nil + El/ iii^CkCew 8\/2iii,, 



^spect; 



;ator 



3/2 



{pkq-piYPi • {pk, -pi) 



xTr 



-1^x75 + )/*(!- 75) (A + m,)75/ 



m 



(86) 



The gamma-matrix factors in this expression that are associated with the B meson corre- 
spond to the leading-twist S- meson light-cone ^bi in Eq. (5). Expanding terms inside the 
trace, we have 

tator -ig^CpCB^l + Ei/m)CKCEW 8^771, 



^specta 



3/2 



{pk,-piypi ■ {pk,-pi) 



X {-Tr y>K^* (1 - 75) M] - (m,/m)Tr (1 - 75) i']} 

_ -ig^CpCBil + Ei/m)CKCEw 8x/2me 



3/2 



{Pkg - PlYPl ■ {Pkg - Pl) 



(Ti + Ts). (87) 
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In evaluating the sizes of the contributions to ^^spectator^ make use of the orders of 
magnitude of the components of the various four vectors in the S-meson rest frame. Some 
of these are given in Eqs. (15a), (18), (25), and (88). In the case of q' , we see from Eq. (22) 
that, in the i?- meson rest frame, 

(g')^ ~ vrric, (88a) 
(g')~ ~ vml/rrib, (88b) 
q'j^ ~ vrric. (88c) 

This is in contrast with either P or q, which, as can be seen from Eq. (25), have plus com- 
ponents that are of order and vrrib, respectively. The suppression of the plus component 
of q' is a consequence of the soft cancellation. 

Now consider the contribution of 71 in Eq. (87). Ti comes from the —'^k^5 term in 
the first parenthesis in Eq. (86), which corresponds to the leading- twist hght-meson light- 
cone distribution. The contribution of Ti is proportional to the one that was considered in 
Ref. [39]. Evaluating the trace in Ti, we obtain 

Ti = -4 {pK ■ e*pb -q'-pK- Pbe" -q'+px- q'e* ■ Pb - ie'''^''Pkae;Pbpq'^) , (89) 

where we have used the convention Tr^y^^jb^^] — Aiea^r^sa'^h^ c' , with 60123 = — e°-^^^ = 
1. It is now easily seen that the terms in Eq. (89) are of order vm\, vrricml, vml, and 
vrricml, respectively. We can compare this contribution with the individual contributions 
that appear before wc apply the soft cancellation or with the contributions in which the 
gluon in Figs. 12(e) and 12(f) attaches at its upper end to the quark lines from the B 
meson or the K meson. In the cancelling contributions, q' in Ti is replaced with P, and Ti 
becomes of order ml/mc- Thus, we see that Eq. (89) is suppressed as vmc/rrib relative to the 
cancelling contributions, in agreement with what we expect for the soft cancellation from 
the general factorization proof. In the contributions in which the gluon in Figs. 12(e) and 
12(f) attaches at its upper end to the quark lines from the B meson or the K meson, which 
are factorizing contributions that contribute to the 5-meson-i^-meson form factor, q' in Ti 
is replaced with pf, or pk- With this replacement, Ti is again of order ml/mc- However, 
as we have mentioned, the Born-level factorizing contributions to the production of a ^Pj 
charmonium with J = or J = 2 vanish, and so the contributions of Eq. (89) are suppressed 
as rricV / {asfrib) relative to the leading factorizing contributions. 
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Next consider the contribution of T2 in Eq. (87), which comes from the 'pKikl^/fn term 
in the first parenthesis in Eq. (86). The leading contributions in T2 are of order vm\ and are 
suppressed by a factor rricV / {asTrib) relative to the factorizing terms. However, the leading 
contributions in T2 are proportional to the transverse components of g^. Upon integration 
of T2 over the angles of the transverse components of Qk, these leading contributions vanish. 
From Eq. (18), we see that the minus component of Qk vanishes and that the plus component 
of Qk is suppressed by a factor Aqcd/?716 relative to the transverse components. Hence, the 
contributions of T2 are suppressed by a factor AqcbttIcV / (cksItiI) relative to the factorizing 
terms, and are negligible in comparison with the other factorization- violating contributions. 

Now let us retain only the leading contribution to ^spectator ^ which is proportional to Ti. 
As wc have mentioned, this contribution is the one that was considered in Ref. [39]. In 
that calculation, light-quark masses were taken to be zero, Qk — ypx — Pkq was neglected in 
comparison with ypK, and pi was taken to have only a plus component, which is written as 
p'l — ^Pq ~ ^p^. Under these assumptions, pk^ — ypK, which has only a minus component 
that is nonzero, q' — qi—Pi{qi-pK)/iPi-PK), and {q')~^ — 0. Then, the resulting contribution 
is 

^spectator ^ -^^^^^^^(1 + i^^/m)C^CEW l^V^m, 

nI'"^ iy'^PK ■ PhPi ■ pk 

X {pk ■ e*Pb -q' -Pk- Pte* ■ q' - ie'"'^^PkaelPb0'^ , (90) 

which yields an endpoint divergence, owing to the factor y"^ in the denominator. The terms in 
parentheses in Eq. (90) are all of order vrricml. That is, there is a suppression factor, relative 
to the cancelling terms, of order ml/ml. From the arguments of Sec. HIK, we expect such 
a suppression because we are neglecting the transverse momenta of the constituents of the 
B meson and the hght meson. In order to make contact with the calculation of Ref. [39], we 
make the replacements CBi^ + Ei/m) — ;> {—i/A)fB^B{0 and Ck — > {i / ^) f K{y) ■, multiply 
by a factor A/2m^^/(2mc)2 ^ ^/TJm'c to compensate for the normalization of the quarkonium 
state relative to the normalizations of the quark and antiquark states, multiply by a P-wave 
quarkonium spatial wave function, and integrate over the wave-function momentum. Then, 
decomposing q and e* into J = and J — 2 angular-momentum tensors, we obtain agreement 
between Eq. (90) and the divergent terms in Eqs. (15) and (17) of Ref. [39].^^ We find in 

We have also checked that, if we keep the exact expression, rather than taking the soft approximation, 
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the J — 1 case that 

= -^w^ L — /. * '''' 



where /n is defined in Ref. [39] and we have retained only the infrared-divergent terms. 



V. SUMMARY AND DISCUSSION 



In this paper, we have given detailed proofs, valid to all orders in as-, of factorization the- 
orems for two exclusive quarkonium-production processes: the production of two quarkonia 
in e+e~ annihilation and the production of a charmonium and a light-meson in S-meson 
decays. We have supplemented our proofs with one-loop examples of the factorization and 
cancellation of soft singularities. (See Sec. IV.) Proofs of these factorization theorems were 
sketched in Ref. [15]. In the present paper, we have provided more detailed arguments. 
The proofs in Ref. [15] did not consider the possibility that on-shell lines could emit gluons 
with arbitrarily small momenta. Such a possibility arises, for example, when one computes 
short-distance coefficients by making use of on-shell matching conditions. In the present 
paper, we have shown that factorization still holds when one takes into account this possi- 
bility. We have also given more refined estimates of the violations of factorization than were 
given in Ref. [15], by considering the dependence of such violations on the velocity v of the 
heavy quark or antiquark in the quarkonium rest frame. We note that, although our proofs 
are demonstrated in models in which external lines are taken to be on the mass shell, the 
methods of these proofs would apply to off-shell models as well, provided that the models 
maintain gauge invariance. 

In the proofs of factorization, our general strategy has been to identify soft singularities, 
collinear singularities, and would-be collinear singularities that appear in the limit of zero 
heavy-quark mass. By demonstrating the factorization of these singularities and would-be 
singularities, we are able to argue that the associated logarithmic enhancements also factor- 
ize. Once the logarithmic enhancements have been removed, the remainder of the production 
amplitude can depend only on the hard scale and, hence, is perturbatively calculable. 

In demonstrating the factorization of singularities and would-be singularities, we have 
made use of standard techniques (see, for example, Rcfs. [25-28]), but we have had to 

then we obtain the finite terms in Eqs. (15) and (17) of Ref. [39]. 
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augment them in order to deal with the situation in which low-energy coUinear gluons 
attach to soft gluons. For this purpose, we made use of the approach developed in Ref. [16] 
in the context of the production of light mesons in e'^e~ annihilation. The methods of proof 
that we have given here should, generally, be applicable to proofs of factorization for other 
exclusive processes in QCD. 

Our factorized form for exclusive production of two quarkonia in e+e" annihilation is given 
in Eq. (1). The expression in Eq. (1) has been used in leading-order and next-to-leading- 
order calculations of exclusive double-charmonium production. It is generally referred to as 
the "NRQCD factorization" formula. 

Our factorized form for the exclusive production of a charmonium and a light meson 
in B decays is given in Eq. (2). An expression of the form in Eq. (2) was suggested in 
Ref. [12] on the basis of an analysis of B-meson decays to light mesons. (In Ref. [12], the 
factorized form was written in terms of the quarkonium light-cone distribution, rather than 
in terms of NRQCD matrix elements.) In Ref. [12], it was conjectured that the violations 
of factorization should vanish in the limit rric — > 0, but a detailed analysis of the scaling of 
the violations of factorization with rric, mb, and v was not given. 

We find, generally, that the violations of factorization arc suppressed by a factor /j = 
rricv'^/Q for each charmonium i in an 5"- wave state and by a factor /j = rndQ for each 
charmonium i in a higher orbital-angular-momentum state, where Q — y/s is the CM energy 
in e'^e~ annihilation and Q = in S- meson decays. Because the violations of factorization 
are proportional to v, they vanish (up to corrections that are proportional to Aqcd/<5) if one 
works to order v'^ in one charmonium. This statement has been confirmed in calculations at 
order (Refs. [13, 37]). 

In the case of S-meson decays, the error-suppression factors for S'-wave and P-wave 
charmonia are fi — v'^mc/rri}, and /j = rric/mb, respectively. These are not particularly small, 
and the violations of factorization may well be comparable to the factorized contributions. 
In the case of e+e" annihilation, the error-suppression factors are smaller by a factor nib/y/s 
than in the case of S-meson decays. Furthermore, there is a suppression factor for each 
quarkonium in the process. Hence, in the case of e+e~ annihilation, the errors are likely to 
be sufficiently small that the factorization formula would be useful. Since the coefficients of 
the suppression factors are nonperturbative quantities, their sizes must be determined, at 
present, through phenomenological studies. 
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In special cases, the relative sizes of the violations of factorization may be enhanced be- 
cause of quantum-number considerations. For example, in S-meson decays, the production 
of ^Pq or charmonia through factorized contributions is not allowed in order a^. The 
production of ^Pq or ^P2 charmonia through factorization- violating contributions does occur 
in order a^. Therefore, in these cases, the violations of factorization are enhanced by a 
factor l/cts relative to the factorized contributions. 

Finally, we mention that we could have written the coUinear functions J"^ that are asso- 
ciated with each quarkonium in terms of light-cone distributions instead of NRQCD matrix 
elements. As we have explained in Sec. IIIK, such an approach yields a hard-scattering 
function H that is manifestly free of soft divergences. In contrast, in the factorized ex- 
pressions involving NRQCD matrix elements that we have presented, power-suppressed soft 
divergences do appear in H and must be discarded. In the case of double-charmonium 
production, these soft divergences are suppressed as fifj, while the corrections to the light- 
cone-distribution factorization formula are suppressed only as fi + fj. Therefore, in the 
case of double-charmonium production, the factorized expression involving NRQCD matrix 
elements that we have presented is more accurate than a factorized expression involving 
light- cone distributions for the quarkonia. 
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Appendix A: Spin projectors 

In this appendix we derive quark-antiquark spin projectors for the case in which the 
quark and antiquark have different masses. We take the momentum of the quark to be Pg 
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and the momentum of the antiquark to be Pg, with both the quark and the antiquark on 
shell: — nig, — m?. Therefore, in the quark- antiquark CM frame we have 

Pg^{Eg,+q), (Ala) 
Pg^{Eg,-q), (Alb) 

where Eg = ^vrig + g^, Eg = + q^, and q is the quark three-momentum in the CM 

frame. It follows that 

P, = + «• (A2a) 

where P — Pg + Pq and, in the CM frame, P and g are given by 

P = (M,0), (A3a) 
g=(0,q), (A3b) 

where M = Eg + Eg. 

The quark and antiquark spinors are given by 

u{Pq, Sg) = Ng\ I , (A4a) 



-q ■ (Tr]{sg] 
{Eg + mg)rj{sgj 



V(Pg,Sg)^Ng\ ^'"^^'^"^ I, (A4b) 



where the normalization factors are 



, \2Ei(Ei + rrii)] ^nonrelativistic, 
A^i = r ^ /^^ ' , (A5) 

[Ei + rrii] ^relativistic, 

for i = g or q. Here, ^(sg) and r]{sg) are the two-component spinors for the spin states Sg 
and Sg, respectively, with 77 in a representation that is conjugate to that of ^. It then follows 
straightforwardly that the spin-singlet projector is given by 

^l{Pq,mg,Pg,mg) = ^ U{Pg , S g)V {Pg , S g) {^S g, g\00) 

N N- 

' ' :{i>, + rng){r+Eg + Eg)^^{jt>g-mg) 



2{Eg + Eg)V2 

-I- 'm- \ M M- 

(A + mg)75(/^g -m^), (A6) 



{Eg +mg + Eg + Ulg) NgNg 



2V2{Eg + Eg) 
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and the spin-triplet projector is given by 



^3{Pq,mq,Pq,mq,X) = ^ u{pq, Sq)v{pq, Sq) (ISq, ^SqllX) 



{i>q + rUq) (f + Eq + Eq) / (A) - TUq) , (A7) 



2{Eq + Eq)V2 

where e(A) is the polarization vector of the spin-triplet state whose components in the 
quarkonium rest frame are 



€(±) =T^(l,±i,0), (A8a) 
e(0) = (0,0,1). (A8b) 

The results in Eqs. (A6) and (A7) are equivalent, in the equal-mass case to those in Ref. [40]. 

Note that the spin projectors Ilj in Eqs. (A6) and (A7) are for the decay of a qq pair. 
The projectors Ilj for the production of a qq pair can be obtained in a similar manner as 



^liPq,'mq,Pq,mq) = ^ (| Sg, ^ Sq\00)v{pq, Sq)u{pq, Sq) 



N N- 

' ' {i)q-mq)^^{r+Eq + Eq){i)q + mq) 



2{Eq + Eq)V2 
(Eq + mq + Eq + TTlq) NqNq 



2^f2{Eq + E-^ 

^3(Pq,mq,Pq,mq,X) = ^{^Sq, ^ Sq\lX)v(pq, Sq)u(pq, Sg) 



{i)q-mq)-i^{i)q + mq), (A9a) 



N N- 

' ' {^^-mq)^*{X){r+Eq + Eq){^q + mq). (A9b) 



2{Eq + Eq)V2 

The relationship between Ilj and flj is 

n, = T^n^o. (Aio) 
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